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ment  experiments  was  conducted  using  a  one-fifth  scale  model  of  a  single 
wheelset  on  tangent  track  subjected  to  static  and  dynamic  loading  conditions. 

The  results  of  these  experiments  were  compared  to  simulations  based  on  a 
nonlinear  theory  developed  to  represent  the  important  phenomena  associated 
with  dynamic  wheelclimb. 
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from  quasi-steady  analysis  of  wheel/rail  forces.  Evidence  has  been  found  that 
derailment  criteria  employing  variables  measured  in  addition  to  wheel  loads 
may  be  successful  in  predicting  derailment  safety,  and  that  diagnostic  criteria 
may  be  developed  for  warning  of  impending  derailment. 
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PREFACE  AND  ACKNOWLEDGEMENTS 

In  support  of  the  Improved  Track  Structures  Research  Program  of  the 
Office  of  Research  and  Development  of  the  Federal  Railroad  Administration 
(FRA) ,  the  Transportation  Systems  Center  (TSC)  is  conducting  analytical  and 
experimental  studies  of  the  relationship  between  track  design  and  maintenance 
parameters  and  the  safety  and  performance  of  trains,  resulting  from  the  dynamic 
interaction  of  the  track  and  train.  The  objectives  of  this  effort  are  directed 
towards  the  development  of  a  technological  data  base  that  can  be  used  by 
Government  and  industry  to: 

a.  Establish  safe  operating  regimes  based  on  the  dynamics  of  track/ 
train  interaction. 

b.  Evaluate  the  safety  and  economics  of  various  types  of  train  operation 
and  track  maintenance  practices. 

c.  Formulate  recommendations  for  reducing  derailments  due  to  track/ 
train  interaction. 

The  approach  being  employed  in  the  conduct  of  these  studies  and  in 
development  of  engineering  guidelines  to  meet  the  above  objectives  includes 
the  following: 

a.  Definition  of  performance  indices  which  are  related  to  the  likelihood 
of  vehicle  derailment,  lading  damage,  passenger  ride  quality,  and 
operational  efficiency. 

b.  Compilation  of  analysis  tools  and  development  of  a  methodology  for 
relating  track,  rail  car,  locomotive,  and  train  parameters  to  the 
performance  indices. 
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c.  Parametric  studies  relating  the  performance  indices  to  track  con¬ 
struction  and  maintenance  parameters  for  the  range  of  in-service 
train  configurations,  equipment  parameters  and  speeds. 

d.  Development  of  techniques  and  methods  for  track/train  dynamics, 
laboratory  and  full-scale  test  design,  and  rail  vehicle  model 
validation,  by  which  the  performance  indices  can  be  quantitatively 
evaluated  and  adjusted  as  required  to  provide  their  proper  threshold 
values . 

The  principal  safety-related  problem  currently  facing  the  railroads  is 
derailment.  Derailments  result  in  direct  loss  due  to  damage  to  rolling  stock, 
track,  lading  and  facilities,  and  in  injury  to  passengers,  and  operators,  and 
liability  for  release  of  hazardous  materials.  Less  well -publicized,  but  of 
great  significance,  is  the  indirect  economic  loss  due  to  decreased  operating 
speeds  and  other  factors  affecting  system  performance  that  are  found  necessary 
because  of  the  probability  of  derailments  occurring.  Advances  in  understanding 
of  derailment  phenomena  would  then  have  the  dual  benefit  of  (1)  improving  the 
technology  to  be  more  resistant  to  derailment,  and  (2)  provide  a  rational 
basis  for  formulating  economically  sound  safety  standards  for  maintenance  and 
operation  that  are  designed  to  reduce  derailment  probability. 

Because  of  the  compelling  importance  of  derailment  as  a  problem  both  to 
the  railroads  and  to  our  national  transportation  system,  and  the  relatively 
undeveloped  state-of-the-art  in  understanding  the  problem,  the  FRA  has  initiated 
a  major  research  program  in  this  area.  This  multifaced  program,  coordinated  by 
TSC,  involves  the  characterization  of  vehicles  and  trackage  previously  described, 
as  well  as  formulation  of  appropriate  models  for  quantifying  derailment  phenomena 
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and  estimating  resultant  derailment  probabilities.  It  is  essential,  however, 
that  data  be  made  available  on  the  important  derailment  phenomena,  so  that 
the  quantitative  models  used  are  soundly  based.  Until  recently,  experiments 
required  to  meet  this  requirement  were  prohibitively  costly,  time-consuming, 
hazardous,  and  potentially  inaccurate  due  to  the  difficulty  in  controlling  test 
conditions.  Recently  techniques  and  apparatus  have  been  developed  to  perform 
meaningful  experiments  on  derailment  through  the  use  of  dynamically  scaled  models. 
In  this  study  many  of  the  requirements  for  quantification  of  derailment  phenomena 
essential  to  the  TSC  program  have  been  met  by  an  experimental  program  which 
utilizes  these  techniques  and  apparatus. 

The  overall  objective  of  this  research  program  is  to  provide  validation 
of  analytical  models  of  wheel/rail  interaction  and  of  performance  indices  that 
can  be  used  to  predict  derailment.  The  experimental  results  from  the  scale 
model  testing  will  be  used  as  follows: 

1.  In  confirming  or  improving  analytical  models  of  the  mechanics  of  . 
wheel/rail  interaction,  especially  under  extreme  conditions  of 
load  or  geometry. 

2.  In  generating  critical  threshold  values  of  performance  indices 
which  will  accurately  signal  impending  derailment. 

3.  In  developing  empirical  models  of  wheel/rail  interaction,  which 
can  be  used  without  recourse  to  the  detailed  theory,  and  which 
may  therefore  provide  a  more  efficient  approach  to  analysis  and 
simulation. 
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During  the  progress  of  this  study  many  useful  technical  comments  and 
suggestions  have  been  contributed  by  staff  at  the  U.S.  Department  of  Transportation, 
Transportation  Systems  Center,  Cambridge,  Massachusetts.  Particularly  significant 
*  were  the  contributions  of  Mr.  H.S.H.  Lee,  the  Technical  Monitor  of  the  contract. 

Dr.  Russel  Brantman,  and  Mr.  William  Thompson,  III.  The  authors  acknowledge 
the  contributions  of  Princeton  University  technical  staff  W.F.  Putman,  E.  Griffith, 

®  and  J.  Grieb,  and  undergraduate  students  J.  Karohl  and  P.  Moy.  The  report 
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EXECUTIVE  SUMMARY 


Introduction 

Statement  of  Problem: 

The  most  widely  accepted  criterion  for  wheelclimb  derailment  defines  an 
upper  limit  on  wheel/rail  contact  forces  on  the  climbing  wheel,  the  limit 
varying  with  the  time  duration  of  the  forces.  The  reliability  of  this  criterion 
in  predicting  derailment  is  of  major  importance,  since  it  is  used  to  evaluate 
the  derailment  safety  of  railroad  vehicles,  operating  conditions,  and  track 
structures.  This  report  presents  an  analytical  and  experimental  evaluation  of 
time-dependent  wheel  load  derailment  criteria  for  wheelclimb,  which  shows  that 
this  type  of  criterion  alone  is  inadequate  for  derailment  prediction. 

Main  Objectives  of  Current  Research: 

The  central  goal  of  this  program  was  to  develop  a  fundamental  understand¬ 
ing  of  the  wheelclimb  derailment  process.  Specific  objectives  undertaken  to 
meet  this  goal  were: 

a.  Establishment  of  an  experimental  data  base  of  derailment  events, 
using  a  dynamically  scaled  model  of  a  single  wheelset,  for  use  in 
evaluating  the  validity  of  derailment  criteria.  Evaluation  of 
specific  candidate  time-duration  dependent  wheel  load  criteria, 
such  as  that  proposed  by  the  Japanese  National  Railways  (JNR) . 

b.  Formulation  of  analytical  models  for  predicting  wheelclimb  under 
quasisteady  and  dynamic  loading  conditions. 

c.  Experimental  validation  of  the  analytical  models  using  scale  model 
experimental  data. 
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Technical  Approach: 

Derailment  may  occur  as  a  result  of  several  distinct  processes,  including 
wheelclimb,  wheel  lift,  rail  rollover,  gauge  spreading,  and  component  failure, 
each  requiring  specification  of  allowable  wheel/rail  forces  or  other  measures. 
This  report  concerns  the  first  of  these  processes,  which  is  directly  related  to 
the  dynamics  of  the  vehicle  on  curved  and  tangent  track,  and  to  the  track  mis¬ 
alignments  that  excite  the  vehicle.  Wheelclimb  is  a  derailment  process  in 
which  large  lateral  forces  acting  on  the  wheelset  cause  on  wheel  to  climb  up 
and  over  the  rail.  When  this  wheel  is  in  flange  contact,  large  restoring 
forces  resisting  wheelclimb  result  from  the  large  contact  angle  at  the  wheel/ 
rail  interface.  Frictional  effects,  known  as  creep  forces,  aid  or  inhibit  the 
process,  depending  on  the  angle  of  attack  of  the  wheel  relative  to  the  rail. 

The  vector  sum  of  the  contact  forces  may  be  measured  by  an  instrumented 
wheelset.  Since  the  maximum  allowable  lateral  force  generally  scales  with 
vertical  load  on  the  derailing  wheel,  a  derailment  quotient,  or  L/V  ratio,  is 
often  used  as  a  measure  of  safety  (L  and  V  are  the  lateral  and  vertical  force 
components,  respectively).  Under  dynamic  conditions,  significantly  larger 
derailment  quotients  may  occur  without  derailment  than  would  be  expected  under 
quasisteady  operating  conditions.  The  JNR  derailment  criterion  was  developed 
to  account  for  this  time-duration  dependency. 

To  study  the  wheelclimb  process  in  detail,  three  distinct  processes  were 
identified  (Table  EXEC.l): 

a.  Quasisteady  wheelclimb,  in  which  lateral  velocity  is  negligible  and 
the  yaw  angle  remains  essentially  constant.  This  process  applies  to 
detailment  during  steady  state  curving. 

b.  Single  degree-of-freedom  (DOF)  wheelclimb,  in  which  lateral  velocity 
effects  are  important,  but  the  yaw  angle  remains  essentially  constant. 
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This  condition  describes  steady  state  curving  plus  train  action, 
discrete  changes  in  track  geometry,  or  truck  hunting  with  relatively 
stiff  primary  yaw  suspensions  for  the  wheelset. 

c.  Two  degree-of-freedom  wheelclimb,  which  includes  lateral  velocity 
and  changing  yaw  angle  effects.  This  process  describes  dynamic 
curving,  wheelset  hunting,  and  truck  hunting  with  relatively  compliant 
yaw  primary  suspensions. 

In  this  study  the  three  derailment  processes  have  been  analyzed  theoreti¬ 
cally  in  detail.  Experiments  using  a  scale  model  wheelset  on  tangent  track 
was  used  to  validate  the  theoretical  analysis.  A  computer  program  was  developed 
for  simulation  of  the  wheelset,  including  coupling  to  the  truck  frame,  which 
permits  it  to  be  integrated  readily  into  complete  vehicle  simulations. 

A  comprehensive  series  of  experiments  was  performed  to  evaluate  existing 
and  proposed  derailment  criteria  for  quasisteady  and  single  DOF  dynamic  wheel - 
climb,  and  to  validate  theoretical  models  for  all  three  processes. 

Conclusions 

The  experimental  and  analytical  studies  of  wheelclimb  derailment  presented 
in  this  report  result  in  the  following  conclusions  regarding  derailment  cri¬ 
teria: 

1)  Wheelclimb  derailment  criteria  based  on  quasisteady  theory  are  adequate 
for  derailment  prediction  under  conditions  of  negligible  lateral  inci¬ 
dent  velocity  and  constant  positive  yaw  angle.  At  negative  yaw  angles 
derailment  occurs  at  L/V  ratios  somewhat  below  the  theoretical  limits 
(Figure  EXEC.l). 

2)  Application  of  the  JNR  and  other  time-duration  dependent  derailment 
driterion  for  nonderailment,  marginal  derailment,  and  complete 
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derailment  cases  are  each  spread  over  an  order  of  magnitude  in  L/V 
ratio  amplitude  and  time  duration.  None  of  the  criteria  tested 
could  distinguish  between  safe  and  unsafe  conditions.  Defining  a 
safety  criterion  below  the  derailment  data  would  be  overly  conser¬ 
vative,  possibly  excluding  vehicles,  operating  conditions,  and  track 
maintenance  standards  that  could  be  demonstrated  to  be  safe  with  a 
more  discriminating  criterion.  Furthermore,  the  data  provided  do 
not  necessarily  establish  a  lower  bound  for  derailment,  since  lower 
points  possibly  could  be  measured  (Figure  EXEC. 2). 

3)  The  analytical  models  for  dynamic  wheelclimb  yield  accurate  predic¬ 
tions  of  wheelset  response  to  external  force  inputs,  in  term,  of 
wheelset  motions  during  derailment  and  wheel/rail  interaction  forces 
(i.e.  L/V  ratios)  (Figure  EXEC. 3). 

4)  Evidence  has  been  found  that  derailment  criteria  employing  variables 
measured  in  addition  to  wheel  loads  may  be  successful  in  predicting 
derailment  safety,  and  that  diagnostic  criteria  may  be  developed  for 
warning  of  impending  derailment  (Figure  EXEC. 4) . 

Recommendation  for  Future  Research 

To  achieve  the  objective  of  defining  a  reliable  measure  of  derailment 
safety,  the  following  would  be  useful  as  topics  of  future  research: 

a)  New  Wheelclimb  derailment  criteria  should  be  developed  and  validated 
that  include  variables  in  addition  to  wheel  loads,  such  as  lateral 
velocity  and  yaw  angle,  that  are  readily  measured  under  full  scale 
test  conditions.  Such  multivariable  criteria  should  be  better 
indicators  of  derailment  safety. 
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b)  Criteria  for  the  wheelclimb  process  should  be  combined  with  criteria 
for  other  derailments  modes,  such  as  gauge  spreading  and  rail  rollover, 
to  yield  a  comprehensive  safety  measure.  Track  stiffness  parameters 
would  be  key  variables  in  the  comprehensive  criteria  formulated. 

c)  The  results  of  analytical  and  experimental  studies  of  derailment  of 
single  wheelsets  should  be  extended  to  complete  trucks  and  vehicles. 

It  is  very  important  to  establish  the  degree  to  which  single  wheelset 
criteria  may  be  applied  directly  to  complete  vehicle  configurations. 

In  addition  to  the  above,  the  concept  developed  in  this  study  of  derail¬ 
ment  diagnostics  should  be  pursued.  This  technique  provides  a  means  for 
detecting  marginally  safe  conditions  in  full  scale  tests,  so  that  safety- 
related  phenomena  may  be  measured  without  actual  derailments  being  required 
in  the  test  plan.  Scale  model  experiments  will  continue  to  be  useful  to  the 
study  of  fundamental  derailment  processes,  due  to  the  greater  control  of  test 
conditions  possible  and  the  relative  ease  of  study  of  the  full  range  of 
derailment  conditions. 
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TABLE 

EXEC.i.  COMPARISON  OF  WHEELCLIMB 

DERAILMENT  PROCESSES 

PROCESS 

IDENTIFICATION 

NAME 

LATERAL 

VELOCITY 

YAK- 

ANGLE 

DESCRIPTION 

A 

Quasisteady 

wheelclimb 

Negligible 

Constant 

Steady  state  curving 

B 

i 

Single  DOF 
wheelclimb 

Important 

Constant 

Steady  State 
curving  plus  train 
action,  discrete  changes 
in  track  geometry,  truck 
hunting  with  stiff  primary 
yaw  suspensions 

C 

Two  DOF 
wheelclimb 

Important 

Changing 

Dynamic  curving,  wheel  set 
hunting,  truck  hunting  with 
compliant  yaw  suspensions 

L 


J 


JNR-CRITERION 

NON-DERflILMENTS 

NOTATION  FOR  TAN  ANGLEi 


X>  >3  DEG 
Bi  0  DEG 
♦.  3  OEG 


JNR  CRITERION 

• 

i 

\  . 

8  e  • 

♦  ®  s 

go* 

©  J? 

+ 

> 

X 

+  + 

T1  (SEC) 


Figure  EXEC. 2a  Comparison  of  Dynamic  Wheelclimb  Experimental  Data  with 
JNR  Criterion  for  Nonderailment  Cases. 
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Figure  EXEC. 2b  Comparison  of  Dynamic  Nheelcliab  Experimental  Data  with 
JNR  Criterion  for  Marginal  n*railr»«t  C«se«. 
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1 .  INTRODUCTION 

The  most  widely  accepted  criterion  for  wheelclimb  derailment  defines  an 
upper  limit  on  wheel/rail  contact  forces  on  the  climbing  wheel,  the  limit 
varying  with  time  duration  of  the  forces.  The  reliability  of  this  criterion 
in  predicting  derailment  is  of  major  importance,  since  it  is  used  to  evaluate 
the  derailment  safety  of  railroad  vehicles,  operating  conditions,  and  track 
structures.  This  report  presents  an  analytical  and  experimental  evaluation  of 
time -dependent  wheel  load  derailment  criteria  for  wheelclimb,  which  shows  that 
this  type  of  criterion  alone  is  inadequate  for  derailment  prediction. 

To  study  wheelclimb  derailment  processes  and  evaluate  derailment  criteria, 
a  series  of  derailment  experiments  was  conducted  using  a  one-fifth  scale  model 
of  a  single  wheelset  on  tangent  track  subjected  to  static  and  dynamic  loading 
conditions.  The  results  of  these  experiments  were  compared  to  simulations 
based  on  a  nonlinear  theory  developed  to  represent  the  important  phenomena 
associated  with  dynamic  wheelclimb. 

This  study  shows  that  the  JNR  and  other  time-duration  dependent  criteria 
based  on  wheel  load  measurements  alone  are  unsuccessful  in  predicting  derailment 
safety  for  dynamic  wheelclimb,  For  wheelclimb  processes  involving  negligible 
lateral  velocities,  the  derailment  limit  can  be  estimated  from  quasisteadv 
analysis  of  wheel/rail  forces,  Evidence  has  been  found  that  derailment  criteria 
employing  variables  measured  in  addition  to  wheel  loads  may  be  successful  in 
predicting  derailment  safety,  and  that  diagnostic  criteria  may  be  developed  for 
warning  of  impending  derailment. 
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There  is  a  strong  aotivation  for  using  wheel  load  measurements  for  de¬ 
railment  prediction.  The  mechanics  of  derailment  for  a  complete  truck  are 
complex  and  highly  nonlinear.  At  present  no  analytical  model  for  complete 
trucks  has  been  developed  and  validated  that  can  relate  measured  variables 
other  than  wheel  loads  to  derailment  occurrence.  The  hope  that  wheel  loads 
would  be  useful  to  derailment  research  has  led  to  the  development  of  instrumented 
wheelsets  as  major  experimental  research  tools,  which  are  also  employed  to 
study  problems  associated  with  failure,  wear,  and  large  deformation  of  wheels 
and  track  [22-24].* 

The  need  for  accurate  and  reliable  methods  for  derailment  prediction  has 
intensified  recently  as  efforts  increase  to  relate  vehicle  characteristics, 
train  operation,  and  track  maintenance  standards  to  system  safety.  Kith 
increasing  success  the  relationships  between  vehicle  and  track  descriptors  and 
the  resulting  wheel/rail  reaction  forces  have  been  quantified  using  analytical 
or  empirical  means  [24-28] .  The  critical  missing  links  are  definitive  specifications 
relating  these  forces  to  safety. 

Derailment  may  occur  as  a  result  of  several  distinct  processes,  including 
wheelclimb,  wheel  lift,  rail  rollover,  gauge  spreading ,  and  component  failure, 
each  requiring  specification  of  allowable  wheel/rail  forces  or  other  measures  [28]. 
This  report  concerns  the  first  of  these  processes,  which  is  directly  related  to 
the  dynamics  of  the  vehicle  on  curved  and  tangent  track,  and  to  the  track  mis¬ 
alignments  that  excite  the  vehicle.  Wheelclimb  is  a  derailment  process  in  which 
large  lateral  forces  acting  on  the  wheelset  cause  one  wheel  to  climb  up  and  over 
the  rail  (see  Figure  1.1).  When  this  wheel  is  in  flange  contact,  large  restoring 
forces  resisting  wheelclimb  result  from  the  large  contact  angle  at  the  wheel/ 
rail  interface.  Frictional  effects,  known  as  creep  forces,  aid  or  inhibit  the 

•Numbers  in  brackets  refer  to  the  list  of  References. 
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process  depending  on  the  angle  of  attack  of  the  wheel  relative  to  the  rail. 

The  vector  sum  of  the  contact  forces  are  measured  by  an  instrumented  wheelset, 
and  are  distinct  from  the  forces  applied  to  the  wheelset  through  its  bearings. 

Wheelclimb  may  occur  at  low  forward  and  lateral  velocities,  with  the 

resultant  wheel/rail  forces  predicted  using  a  quasisteady  approach  [5,8]  ,  and 

have  been  confirmed  experimentally  in  [8,30,31].  Since  the  limit  on  lateral 

forces  F  scales  with  vertical  axle  load  F  a  derailment  quotient,  or  L/V 
^o  zo 

ratio,  is  often  used,  where  L  and  V  are  the  lateral  and  vertical  components  of 
the  wheel/rail  reaction  forces.  When  significant  lateral  velocities  occur  during 
the  wheelclimb  process,  lateral  forces  exceeding  the  quasisteady  derailment  limit 
generally  occur,  but  only  for  relatively  short  time  durations  [24-26].  A 
criterion  developed  by  the  Japanese  National  Railways  (JNR)  attempts  to  account 
for  the  large  contact  forces  by  increasing  the  limit  on  allowed  L/V  ratios 
proportionally  with  1/tj,  where  tj  is  the  duration  of  the  force  pulse  when 
it  is  less  than  50  msec.  [32]. 

The  objective  of  the  analytical  and  experimental  research  summarized  in  this 
report  is  the  evaluation  of  the  validity  of  wheelclimb  criteria  based  on  wheel 
load  measurements.  In  order  to  systematically  treat  this  objective,  three 
distinct  wheelclimb  derailment  processes  may  be  identified,  which  become  the 
basis  for  detailed  analysis: 

■  Process  A  -  Quasisteady  wheelclimb,  in  which  lateral  velocity  is 
negligible  and  the  yaw  angle  remains  essentially  constant.  This 
process  applies  to  derailment  during  steady  state  curving. 

•  Process  B  -  Single  degree -of- freedom  wheelclimb,  in  which  lateral  velocity 
effects  are  important,  but  yaw  angle  remains  essentially  constant. 

This  condition  describes  steady  state  curving  plus  train  action, 
discrete  changes  in  track  geometry,  or  truck  hunting  with  relatively 
stiff  primary  yaw  suspensions  for  the  wheelset. 
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Process  C  -  Two  degree-of-freedom  wheel  climb,  which  includes 
lateral  velocity  and  changing  yaw  angle  effecis.  This  process 
describes  dynamic  curving,  wheelset  hunting,  and  truck  hunting 
with  relatively  compliant  yaw  primary  suspensions. 

In  this  study  the  three  derailment  processes  described  above  are  analyzed 
theoretically  in  detail.  Experiments  using  a  scale  model  wheelset  on  tangent 
track  are  used  to  validate  the  theoretical  analysis.  A  computer  program  is 
developed  for  simulation  of  the  wheelset,  including  coupling  to  the  truck 
frame,  which  permits  it  to  be  integrated  readily  into  complete  vehicle  simu¬ 
lations.  A  comprehensive  series  of  experiments  is  described  which  are  used 
to  evaluate  existing  and  proposed  derailment  criteria  for  quasisteady  and 
single  degree  of  freedom  dynamic  wheelclimb  (Processes  A  and  B) . 

Section  1  describes  the  motivation  for  this  study,  and  reviews  present 
knowledge  in  the  field.  Section  2  presents  the  analytical  development  of 
criteria  for  quasisteady  wheelclimb.  The  experimental  data  presented  agrees 
well  with  the  derived  criterion.  In  Section  3,  a  nonlinear  model  is  validated 
for  wheelset  displacement  responses  to  pulses  in  applied  lateral  forces. 
Simulations  using  the  dynamic  model  indicate  that  the  magnitude  and  duration 
of  L/V  ratios  that  would  be  measured  during  dynamic  wheelclimb  are  not 
related  in  a  functional  way  to  derailment  occurence. 

In  Section  4  time-duration  dependent  wheel  load  criteria  for  dynamic 
wheelclimb  are  evaluated.  A  variety  of  wheel  load  criteria,  including  the 
widely  recognized  JNR  criterion,  are  applied  to  wheel  load  measurements  taken 
from  112  derailment  events,  with  none  found  to  be  successful  in  predicting 
derailment  safety. 
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Appendix  A  summarizes  the  results  of  application  of  the  variety  of  time-  j 

duration  dependent  wheel  load  criteria  to  the  112  derailment  events  recorded  j 

during  the  experiments.  Appendix  B  describes  in  detail  the  apparatus  used  j 

i 

during  Phase  One  of  the  research  program,  with  the  Phase  Two  apparatus  described  j 

in  Appendix  C.  The  computer  simulation  program  used  to  model  wheelclimb 
derailment  is  listed  in  Appendix  D. 
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2.  QUASISTEADY  WHEELCLIMB 

2.1  Analytical  Development 

The  steady  rolling  of  a  wheelset  under  load  is  a  quasistatic  process  in 
which  an  equilibrium  condition  exists  foT  all  applied  forces  and  moments  due 
to  vehicle  and  track  loading.  If  for  a  given  set  of  applied  axle  loads  and 
wheelset  yaw  angle  an  equilibrium  condition  exists,  characterized  by  a 
certain  wheelset  lateral  position,  wheelclimb  derailment  will  not  occur. 

If  no  stable  equilibrium  exists,  wheelclimb  derailment  results.  Yaw  angle 
is  defined  to  be  positive  when  the  derailing  wheel  is  steered  into  the  rail. 

For  convenience,  derailment  is  assumed  to  occur  on  the  right  rail,  viewed 
in  the  direction  of  travel. 

Since  wheelset  yaw  angles  are  generally  small,  force  and  moment  equilib¬ 
rium  conditions  may  be  applied  in  a  vertical  plane  passing  through  the  axle. 

In  this  report  the  effect  of  translation  of  the  wheel/rail  contact  points 
longitudinally  along  the  rails  is  not  considered.  Analysis  of  the  wheelset 
equilibrium  conditions  requires  calculation  of  the  forces  at  the  contact 
points  due  to  longitudinal,  lateral  and  spin  creep,  each  of  which  varies 
with  wheelset  lateral  displacement,  yaw  angle,  and  axle  angular  velocity. 

The  procedure  for  calculating  the  highly  nonlinear  wheel/rail  contact  forces 
in  flange  contact  is  presented  in  detail  in  [6],  and  is  summarized  schematically 
in  Figure  2.1.  The  approach  is  similar  to  that  presented  in  [5]  and  [13], 
each  analysis  including  the  Kalker  creep  force  theory  and  iterative  solutions 
for  axle  velocity  and  normal  forces  at  the  wheel/rail  contact  points.  Numerical 
results  are  produced  for  specific  wheel /rail  contact  geometries,  computed  for 
given  wheel  and  rail  profiles  using  algorithms  developed  in  [17].  Creep 
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forces  at  each  contact  point  are  determined  through  interpolation  of 
numerical  results  from  Kalker's  Simplified  Theory,  tabulated  for 
appropriate  combinations  of  nondimen sional  creepages  and  contact  ellipse 
geometry  [14]. 

The  quasisteady  derailment  limit  (Process  A)  is  derived  by  calculating 

the  maximum  applied  lateral  force  that  may  be  sustained  for  a  given 

vertical  force  and  roll  moment  applied  to  the  wheelset.  The  net  wheel/ 

rail  reaction  forces  and  moments  acting  on  the  wheelset  in  steady  state 

as  shown  in  Figure  2.1  are  balanced  by  applied  forces  F  ,  F  ,  and 

yo  zo 

moment  M. 

Lateral  force: 


0  =  F 


*  l  F  (uf  cosa  ♦  sino.) 
yo  i=l  1  1 


(2.1) 


Vertical  force: 

2 

0  *  F  ♦  1  F,  .  (-yf?  .  sino.  ♦  coso.) 
Z0  iml  -5*1  1  1 


Roll  Moment: 


0  « 

*o 


2 

I  (-1)1L.F,  .  (-y  f.  .sino.  ♦  cosa.) 
1  2  0,2  +  9  1  1  1 


(2.2) 


(2.3) 


The  normal  forces  Fj  .  at  the  contact  points  may  be  eliminated  to  yield  the 
vertical  and  lateral  forces  acting  on  each  wheel.  The  wheel/rail  reaction 
forces  are  then  combined  to  give  l/V  ratios  for  the  wheelset  axle  (both 
wheels)  or  for  the  climbing  wheel  only. 
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Parameters  and  T 2  are  the  individual  wheel  L/V  ratios  for  the  left  and 
right  wheels,  respectively.  Parameter  n  is  the  dimensionless  roll  moment, 
representing  the  degree  of  asymmetry  in  axle  load  between  the  extremes  of 
wheel  lift  (with  a  practical  range  of  *1).  In  the  above  equations  f- 
and  f2  2  are  the  dimensionless  lateral  creep  coefficients  (including  spin) 
on  the  left  and  right  wheels,  $  the  wheel  set  roll  angle,  and  and  o2  the 
associated  contact  angles,  defined  to  be  positive  counterclockwise.  The 
nomenclature  and  coordinate  system  is  developed  in  [5,6],  For  the 
quasisteady  process,  parameters  and  T2  are  functions  of  lateral  position 
y  and  yaw  angle  for  each  value  of  <J/  a  maximum  for  L/V  exists,  usually 
close  to  the  value  for  y  which  yields  the  maximum  contact  angle  on  the 
climbing  wheel. 
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The  L/V  ratio  for  the  climbing  wheel  is  often  used  itself  as  a 
derailment  criterion,  and  is  defined  using  the  above  parameters  as* 
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climbing  wheel 


tan  (a2«-4>)  ♦  \if2  2 
1-  V*2<2  tan”(a2«-$) 


(2.5) 


For  large  yaw  angles  Eq.  (2.5)  approaches  the  classical  limit  of  Nadal  [8], 
but  for  the  range  of  conditions  -3°  <  t|<  <+3°  the  nondimens ional  creep 
f2  2  may  vary  considerably  with  wheelset  loading,  wheel/rail  profiles,  and 
yaw  angle,  over  the  range  -1  <  f2  2  <  1.  Under  quasisteady  conditions, 
the  creepages  are  functions  only  of  the  wheel  and  rail  profiles,  wheelset 
displacement,  and  yaw  angle  relative  to  the  rail,  so  that  Eqs.  (2.1)  through 
(2.5)  are  independent  of  forward  velocity  V. 

2,2  Method  of  Computation 

The  method  used  for  computing  the  quasisteady  derailment  limit  is 
the  same  as  that  described  in  detail  in  [6],  with  the  exception  of  the 
procedure  used  to  calculate  creep  forces.  In  [6]  creep  forces  are 
approximated  by  calculating  initially  the  forces  due  to  longitudinal  , 
lateral,  and  spin  creepages  separately,  then  added  vectorially  with 
approximate  adhesion  limits  imposed.  Further  study  has  shown  that  at 
negative  yaw  angles  and  large  spin  creepage  (i.e.  flange  contact)  signi- 
cant  errors  accrue  using  this  approximation,  yielding  derailment  limits 
at  large  negative  angles  that  do  not  agree  with  Nadal 's  limit. 

The  new  procedure  uses  creep  force  calculations  bssed  on  Kalker's 
Simplified  Theory,  so  that  creep  forces  in  the  lateral  and  longitudinal 
directions  are  computed  directly  from  the  three  creepage  components  (18]. 
♦For  simplicity  assume  right  wheel  is  climbing. 
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To  make  the  computation  efficient,  nondimen sional  creep  forces  are  pre¬ 
computed  and  stored  using  the  following  procedure: 

1)  The  lateral  displacement  at  which  the  contact  angle  is  maximum 
is  determined.  This  point  is  very  close  to  the  point  of  maximum 
lateral  force  (i.e.  derailment  limit)  and  eliminates  the  need  to 
search  over  the  y-direction. 

2)  The  non-dimensional  contact  ellipse  geometry  factors  a  ■  a/ /ab 
and  F  «  b//ab  are  computed  for  this  lateral  displacement.  These 

factors  are  independent  of  wheel  loads.  The  effect  of  ^  on  these 

parameters  is  neglected,  also  the  nondimensional  spin  creep  x  is  computed. 

3)  For  the  above  values  of  T,  F,  x  and  Poisson's  ratio  the  Kalker 
program  is  run  repeatedly  to  generate  a  table  of  lateral  and 
longitudinal  nondimensional  creep  forces  for  all  possible  values 
of  lateral  and  longitudinal  creepages  and  v^,  respectively. 

The  table  is  expressed  in  coordinates  v  and  a,  which  are  the  polar 
transformations  of  v  and  v  . 

y  * 

4)  For  each  value  of  u  and  v  desired,  the  creep  forces  are 

x  y 

calculated  from  the  table  using  interpolation  based  on  the 
Lagrange  polynomial  method. 

The  resulting  calculations  yield  a  limit  for  large  negative  yaw  angles 
of  L/Vfixle  equal  to  3.1,  which  agrees  with  Nadal's  limit,  in  contrast  to 
a  value  of  2,0  computed  in  [29].  At  large  positive  yaw  angles  both  the 
above  and  previously  used  methods  yield  the  correct  Nadal  limit. 

2.3  Comparison  of  Axle  and  Wheel  L/V  Ratio  Criteria 

Due  to  the  nonlinearity  of  the  creep  force  phenomenon,  wheelclimb 
derailment  limits  based  on  both  axle  and  wheel  forces  vary  with  axle 
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loading.  Using  numerical  results  for  wheels  and  rails  with  new  profiles 
extended  from  [6]  to  cases  with  non-zero  roll  moments,  the  degree  to 
which  each  criterion  may  be  applied  universally  is  demonstrated.  In 
Figure  2.2  wheel  and  axle  L/V  ratios  are  shown  for  varying  axle  vertical 
load.  Both  criteria  are  insensitve  to  vertical  load,  with  minor  variation 
evident  only  in  the  range  -1.4°  <  <  0°.  The  two  criteria  are  related 

for  all  vertical  axle  loads  by  a  single  function  shown  in  Figure  2.3. 

Figure 2.4  shows  the  effect  of  roll  moment  parameter  n  on  the  two  criteria. 

A  positive  roll  moment,  increasing  the  vertical  force  on  the  wheel  in 
flange  contact,  produces  a  larger  axle  L/V  ratio  for  wheelclimb  derail¬ 
ment  for  all  yaw  angles.  The  scaling  of  the  derailment  limit  with  roll 
moment  is  reflected  in  the  relative  insensitivity  in  the  wheel  L/V  ratio. 
Although  the  roll  moment  parameter  n  does  not  appear  explicitly  in 
Eq.  2.5^  it  does  affect  the  solution  for  wheelset  force  equilibrium, 
c?-'Jsing  variations  in  f  .  over  the  range  -1.0°  <  \f/  <  1.0°.  The 

*•  t*- 

relations  between  the  two  criteria  for  varying  roll  moment  is  shown 
in  Figure  2.3. 

Table  2.1  summarizes  the  relative  advantages  of  use  of  criteria  of 
each  type.  For  situations  in  which  axle  loads  including  roll  moments 
are  known,  from  simulation  or  vehicle  measurements,  axle  L/V  ratios  are 
better  since  the  individual  wheel  forces  are  not  required.  If  the  roll 
moment  is  not  known,  as  may  be  the  case  for  field  experiments  on  full  scale 
vehicles,  wheel  L/V  ratios  are  better  when  applied  to  data  from  instrumented 


wheel  sets. 


TABLE  2.1  COMPARISON  OF  AXLE  AND  WHEEL  L/V  RATIO  CRITERIA 


Advantage 

Axle  L/V 

Wheel  L/V 

Insensitive  to  axle  vertical  load 

yes 

yes 

Insensitive  to  axle  roll  moment 

no 

yes 

Shows  variation  with  yaw  angle 

yes 

yes 

Does  not  require  simulation  or  measurement 

of  individual  wheel  contact  forces 

yes 

no 
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2.4  Experimental  Validation  and  Usefulness  of  Model 

The  quasisteady  wheelclimb  experiments*  consist  of  a  lateral  force 
ramp  loading  applied  to  the  wheel set,  with  vertical  axle  load,  yaw  angle 
and  forward  velocity  held  fixed.  The  loading  rate  is  held  small  to 
preserve  quasisteady  conditions.  Typical  time  histories  of  axle  L/V 
ratio  and  lateral  displacement  are  shown  in  Figure  2.5. 

The  time  histories  indicate  a  very  well-defined  derailment  limit. 
Fluctuations  in  the  measured  responses  are  due  to  variations  in  local 
track  geometry;  the  "track  signature"  is  quite  repeatable  when  the 
experiment  is  duplicated  at  the  same  track  location,  even  at  different 
forward  velocities.  At  other  track  locations,  the  fluctuations  over 
the  time  history  differ,  but  the  L/V  derailment  limits  are  consistent 
within  a  range  of  +5%.  The  temporal  frequencies  of  these  fluctuations 
are  sufficiently  low  to  maintain  quasisteady  conditions.  Since  the 
variations  in  track  alignment  and  rail  profile  in  the  scale  model  are 
less  than  that  expected  in  full  scale,  the  fluctuations  in  measured 
responses  shown  are  probably  smaller  than  would  be  expected  in  field 
testing. 

Experiments  such  as  that  illustrated  in  Figure  2.5  were  conducted 
for  the  range  of  yaw  angles  -3°  <_  £  3°.  The  maximum  axle  L/V  ratio 

immediately  before  derailment  for  each  yaw  angle  is  recorded  in  Figure 
2.6  for  the  case  of  zero  applied  roll  moment.  Agreement  between  theory 
and  experiment  is  generally  good,  although  the  theory  overestimates  the 
derailment  limit  for  negative  angles  in  the  transition  region  between 
-0.5°  and  -2.0°. 


*  Apparatus  described  in  Appendix  B. 


Experimental  Tine  Histories  of  Lateral  Displacement  and  Axle  L/V  Ratio  During 
adv  Wheel climb  Derailment. 
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MHEELCLIMB  OEMIULNENT  LIMIT 


TMH  UNCLE.  OECREES 

Figure  2.6  Comparison  of  Quasisteady  Theory  and  Experiment  for  Wheel- 
climb  Derailment  Limit  Kith  Zero  Roll  Moment. 


NHEELCUM8  DERR  ILNENT  LIMIT 
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Figure  2.7  Comparison  of  Quasisteady  Theory  and  Experiment  for 

Wheelclisib  Derailment  Limit  for  Varying  Yaw  Angles  and 
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The  results  for  quasisteady  wheelclimb  are  consistent  with  roller 
rig  data  in  [?]»  where  the  theory  also  overpredicts  the  derailment  limit 
for  negative  yaw  angles  (but  with  a  smaller  error).  The  source  of  the 
discrepancy  in  the  present  tests  has  not  been  positively  identified,  but 
may  be  due  either  to  the  presence  of  the  lateral  track  dynamic  input 
(not  accounted  for  in  theory)  or  to  local  variation  in  friction  coefficient 
or  contact  angle.  Since  the  data  presented  are  the  only  published  results 
for  tangent  track  experiments,  conservative  safety  criteria  should  be 
based  on  the  lower  measured  values  rather  than  the  theory  for  negative 
yaw  angles. 

Three  possible  sources  of  error  may  account  for  this  discrepancy, 
two  experimental  and  one  theoretical.  In  the  vicinity  of  the  derailment 
location  on  the  track,  local  variation  in  rail  profile  curvature  from 
nominal  values  may  produce  significant  changes  in  nondimensional  lateral 
creep  coefficient  f,  Error  in  f.  ,  or  friction  coefficient  p  would 

+  » *  ^9  + 

produce  an  error  in  the  predicted  derailment  limit.  The  third  possible 
source  of  error  results  from  the  high  ratio  of  the  semi-axes  of  the 
contact  ellipse  in  flange  contact.  For  the  wheel  and  rail  profiles  used 
in  the  experiment  the  ratio  a/b  is  16,  which  is  beyond  the  numerical 
range  computed  by  Kalker  and  the  experimental  range  measured  by  Brickie 
IS] .  Therefore  numerical  error  in  the  calculated  creep  forces  is  not 
excluded  as  an  explanation  for  the  difference  between  theory  and  experiment. 

The  effect  of  roll  moment  on  the  measured  axle  L/V  derailment 
limit  is  shown  in  Figure  2.7.  Data  points  shown  without  lines  are  touching 
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the  predicted  derailment  limit  at  the  associated  yaw  angle.  Data  points 
with  lines  indicate  the  differences  between  theory  and  experiment .  The 
application  of  constant  roll  moments  under  rolling  conditions  on  the  track 
was  difficult  experimentally,  resulting  in  some  scatter  in  the  data.  The 
agreement  between  theory  and  experiment  is  good,  with  the  same  differences 
as  discussed  above  for  negative  yaw  angles. 
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3.  DYNAMIC  WHEELCLIMB 


3.1  Analytical  Development 

3.1.1  Introduction 

Numerous  theoretical  models  for  wheelset  dynamics  are  available  in  the 
literature  [1-4].  For  large  amplitude  wheelset  motions,  the  mechanics  and 
dynamics  become  highly  nonlinear,  due  to  the  effects  of  contact  geometry  and 
rolling  friction  creep  forces.  For  motions  during  which  tread  contact  is 
maintained  on  both  wheels,  whether  the  wheel  profiles  are  conical  or  curved, 
the  wheelset  dynamic  equations  include  many  terms  of  approximately  equal 
magnitude.  Under  conditions  of  severe  flange  contact  leading  to  derailment, 
the  relative  magnitudes  of  these  various  dynamic  terms  differ  considerably 
from  their  values  in  the  tread  contact  case,  leading  to  conclusions  different 
in  certain  respects  from  those  of  the  cited  references. 

Compared  to  the  rather  rigorous  analytical  models  of  the  preceeding 
references,  theories  used  for  analysis  of  derailment  mechanics  have  been 
simplified.  Equilibrium  models  for  wheels  and  rails  in  flange  contact  by 
Gilchrist  and  Brickie  [5],  Sweet  [6],  and  Yokose  and  Arai  [7]  have  success¬ 
fully  predicted  quasisteady  derailment  limits.  Dynamic  wheelclimb  models  by 
Sweet  [8]  and  Yokose  [7,9,10]  on  the  other  hand,  have  all  been  based  on 
certain  assumptions  which  have  not  been  fully  verified  by  simulation  or  ex¬ 
periment  . 

In  this  study  a  detailed  model  for  wheelset  derailment  has  been  developed. 
The  wheelset  motion  includes  three  independent  degrees  of  freedom;  lateral 
displacement,  yaw  angle,  and  axle  rolling  angular  velocity.  The  truck 
motion  includes  a  single  degree  of  freedom:  the  lateral  displacement  of  the 


truck  frame. 
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While  nearly  all  analyses  include  the  first  two  notions  of  the  wheelset, 
only  aore  recent  papers  by  Clark  [11]  and  Duffek  [12]  have  accounted  for  the 
deviation  of  the  rolling  velocity  from  the  nominal  value  of  V/rQ  expected  when 
the  wheelset  is  centered  and  rolls  with  zero  slip.  As  shown  in  [6],  during 
steady  flange  contact  the  axle  speed  drops  by  about  3%,  which  should  have  a 
significant  effect  on  calculation  of  the  longitudinal  creepage. 

If  wheel/rail  contact  is  maintained  on  both  wheels,  the  wheelset  verti¬ 
cal  and  roll  degrees  of  freedom  are  not  independent  and  may  be  expressed  as 
a  function  of  the  lateral  motion.  Imposition  of  these  kinematic  constraints 
leads  to  generation  of  significant  dynamic  components  in  the  equations  of 
motion.  A  wheelset  kinematic  model  is  developed  in  this  report  to  account  for 
these  dynamic  effects,  as  well  as  for  computation  of  lateral  creepages  in  the 
contact  zones  as  a  function  of  the  velocity  of  the  wheelset  center  of  mass. 

A  digital  computer  program  for  solving  the  resulting  equations  of  motion 
is  described.  The  program  is  compatible  with  published  programs  used  for 
computation  of  contact  geometry  and  creep  force  functions. 

3.1.2  Kinematics  of  Wheelset  Motion 

Kinematic  Model;  Careful  representation  of  the  kinematics  of  wheelset 
motion  is  needed  for  determination  of  the  creepages  at  the  wheel/rail  inter¬ 
faces  and  the  coupling  of  the  vertical,  lateral  and  roll  motions.  The  kinematic 
analysis  employs  three  sets  of  coordinate  systems,  as  illustrated  in  Figure  3.1. 
The  contact  zone  creepages  depend  on  the  relative  velocities  of  the  contact 
points  of  the  wheels  with  respect  to  the  rails.  The  wheelset  kinematic  model 
is  used  to  determine  these  relative  velocities  in  terms  of  the  velocity  of 
the  wheelset  center  of  mass,  utilizing  the  geometrical  constraints  imposed  by 
the  rails  on  the  wheelset  motion.  The  model  is  also  used  to  compute  the 
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Figure  3.1  -  Definition  of  Coordinate  Systems  for  Wheelset  and  Contact  Zones 
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wheelset  vertical  and  roll  accelerations,  as  functions  of  its  lateral  velocity 
and  acceleration. 

The  wheelset  kinematic  model  is  illustrated  in  Figure  3.2.  A  two-dimen¬ 
sional,  planar  motion  model  is  used  to  derive  the  relationships  between  lateral, 
vertical,  and  roll  motions.  The  instantaneous  center  of  rotation,  point  0,  is 
assumed  to  be  in  the  YZ  plane,  with  its  location  depending  on  the  contact  angles 
and  contact  points.  The  kinematic  relationships  are  derived  by  first  determining 
the  position  of  the  wheelset  center  of  mass,  point  C,  relative  to  the  instantaneous 
center: 


(3.4) 


Figure  3.2 


-  Velocities  Defined  by  Wieelset  Kinematic  Model 


I 


3-6 


♦ 


1 

OC  cot  6 


(3.5) 


Contact  Zone  Velocities 

The  lateral  velocity  in  the  contact  plane  at  each  contact  point  is 
determined  by: 


(Lj*L2 )eoao2  $<*>w  j. 
V2,l  "  ain(|oJ  Yw 

(Lj+LjJcotOj  a«*>w  . 

V2,2  "  ain^lajj+dj)  d)T  *w 


(3.6a) 


(3.6b) 


Creepages  are  defined  using  the  model  of  Gilchrist  and  Brickie  [5],  with 
extensions  to  include  dynamic  effects.  Sign  conventions  have  been  modified 
for  consistency  with  Kalker  theory. 

For  longitudinal  creepage  (non-dimensional): 


-flr  .eos*+V 
V 


♦ 


For  lateral  creepage  (non-dimensional): 


(3.7) 


sin*  4  V2.i 
coa^  V 


For  spin  creepage  (dimensional) : 


3,i 


ama- 


v.i 


(3.8) 


(3.9) 


Similar  expressions  were  developed  by  Elkins  and  Gostling  [13],  although 
they  used  the  wheelset  c.m.  lateral  velocity  Yw  instead  of  V2  The  spin 
creepage  defined  in  Equation  (3.9)  is  a  function  of  contact  geometry  only,  for  a 
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given  value  of  wheel set  lateral  displacement.  The  effects  of  yaw  velocity  and 
slip  on  spin  creepage  have  been  neglected.  For  non-zero  values  of  yaw  angle 
the  wheelset  contact  geometry  and  kinematics  become  three-dimensional,  with 
the  longitudinal  displacements  of  the  contact  points  probably  being  the  most 
important  three-dimensional  effect.  These  effects  have  not  been  included  in  the 
present  analysis. 

Vertical  and  Roll  Accelerations 

The  normal  forces  at  the  wheel/rail  contact  points  depend  on  the  vertical 
and  roll  velocities  and  accelerations  of  the  wheelset.  These  quantities  are 
found  readily  through  application  of  the  chain  rule,  which  relates  them  to  the 
wheelset  lateral  velocity  and  acceleratior.  For  roll. 


d*w  • 

•-**5 


(3.10) 


^  •  2 

w 


♦  -  ttr1  V  .  ♦ 


(3.11) 


Similar  expressions  hold  for  vertical  displacement  Z^. 

The  first  partial  derivatives  are  given  by  Equations  (3.4)  and  (3.5).  Both 
the  first  and  second  derivatives  are  functions  of  contact  geometry  only,  and 
may  be  pre-computed  and  stored  as  functions  of  wheelset  lateral  displacement. 
3.1.3  Creep  Forces  Derived  from  Kalker  Theory 

Geometry  of  Contact  Zone:  The  Kalker  Simplified  Theory  of  rolling  con¬ 
tact  is  used  to  model  the  creepage  forces  at  the  wheel/rail  interface.  It  is 
based  on  the  Hertzian  contact  assumption  [14],  which  expresses  the  contact 
zone  shape  and  dimension  in  terms  of  the  local  geometry  constraints  of  the 
wheel  and  rail.  The  contact  ellipse  semi-axes  a^  and  b^  are  given  by: 
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Normalization  of  Creepages 

Given  the  contact  zone  shape  and  dimension,  and  the  associated  creepages 
and  spin,  Kalker  forms  the  non-dimensional  creep  parameters: 

.  £i_ 

l,i  ycj 

'  ’  ^  ih 


■  “3,1  a 


(3.13) 

(3.14) 

(3.15) 


where  y  is  the  friction  coefficient  and  •  *4^  (the  geometrical  mean  of 
the  contact  ellipse  semi-axes). 

For  any  combination  of  a/b  ratio  and  spin  parameter  X.  Kalker 's  Simplified 
Theory  predicts  the  non-dimensional  longitudinal  and  lateral  creep  forces  f ^  ^ 
and  f2  j.  The  dimensional  creep  forces  applied  by  the  rails  on  the  rolling 
wheel  along  the  contact  zone  system  axes  are: 


Fl,i"  “£l,iV  P3,i  P2,i  *  ”f2,iw  P3,i 


(3.16) 
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In  addition  to  these  forces,  a  pure  moment  is  transmitted,  but  of  negligible 
magnitude. 

The  angular  velocity  of  the  wheelset  (aoving  forward)  will  produce  positive 
spin  on  the  left  contact  zone,  and  negative  on  the  right.  To  account  for  negative 
spin,  the  following  relationships  are  used,  based  on  arguments  of  symmetry  [14], 

*  (3.1 

f  2<  f ,  xWf  2(-l,  •>.  x)-f2(  ,-x)— f2(- 1 

3.1.4  ANALYSIS  OF  WHEELSET  DYNAMICS 
Wheelset  Equations  of  Motion 

The  following  assumptions  are  employed  in  deriving  the  wheelset  equations 
of  motion.  Roll  and  yaw  angles  are  assumed  to  be  small.  Gravity  and  suspension 
forces  are  applied  to  the  wheelset  along  the  XYZ  rail  coordinates  axes.  The 
longitudinal  and  lateral  axes  of  the  contact  ellipse  are  assumed  to  coincide  with 
the  X  axis  and  the  YZ  plane,  respectively.  The  equations  of  motion  are  derived 
using  methods  from  [15]  by  Newtonian  methods  using  Euler  Angles  for  calculation 
of  the  gyroscopic  tdrms. 

The  translational  equations  of  motion  are  derived  in  the  XYZ  rail  coordinate 
system  (see  Figure  3.3).  Forces  and  moments  transmitted  to  the  wheelset  from 
the  truck  frame  through  the  primary  suspension  are  represented  by  vectors  and 
fls,  respectively. 

In  the  Y-direction  (lateral): 


Vw  "£FM<Vf2,iCO,V'inai)  4  FSY 
1 


(3.18) 
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Figure  3.3  -  External  Forces  and  Moments  Applied  to  the  Wheelset 
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In  the  2-direction  (vertical); 


V.  ,inarco,ai>  ♦  Fs 


(3.19) 


The  rotational  equations  of  Motion  are  derived  in  the  rotating  xyz  system. 


About  the  x  axis  (roll): 


53F3,it<"1)i<C°*ai“vf2,iBi,IOi)1,i“(sinVw£2,ieo*ai)rv,i1  4  Max  " 


(3.20) 


About  the  y  axis  (pitch) : 


»2*>. •  vVtvw  * 


(3.21) 


(I  -I  )  *  i  -M 
vtjat  *sx  ay 


About  the  z  axis  (yaw) 


y  /  .  (-l)*f,  .T,  ,L.  "X  i*  ♦  i  *  ♦  (I  -I  -I  )  - 
A-u  1,1  3,1  i  as  rw  VW Ayy  ‘as7 


(3.22) 


Application  of  the  Geometrical  Constraints 

The  geometrical  constraints  of  the  single  point,  continuous  contact 
assumption  can  be  used  to  express  the  normal  forces  at  the  contact  zone  by 
solving  Equations  (3.19)  and  (3.20). 


L 


Denote : 


A  ■  m  Z  -  Fc_ 
w  w  SZ 


B  -  I„+w  ♦  (I„-1yy>  *w  ®  * w  ‘  M.x 

*3,1  "  (Lj+Lj) Ccoeoj-Uf 2 , l*inal ^ 


(3.23) 


(3.24) 


(3.25) 


3,2  *  (Lj+Lj) (cosa2-pf2  2sina2) 
where  f  .  (»Vfl*rV(2)/2 


(3.26) 


Coupling  to  Truck  Frame 

Loads  from  the  truck  frame  are  transmitted  to  the  wheel set  through  the 
primary  suspension  (Figure  3.4).  To  obtain  reasonable  results  for  computed 
normal  forces,  the  truck  frame  lateral  displacement  must  be  included  in  the 
simulation  model.  In  this  study,  the  truck  yaw  and  roll  angles  and  vertical 
displacement  are  treated  as  fixed,  but  these  variables  are  included  in  the 
equations  below  for  generality  and  ease  of  extension  to  future  complete  truck 


models. 


For  a  passenger  car,  the  primary  suspension  is  linear  so  that 


'  yvv  *  yvv 


'«z  •  yvv  ♦  vvv 


(3.2?) 
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VERTICAL 


Figure  3.4  -  Model  used  for  single  wheelset,  primary  suspension, 
and  truck  frame  motion. 
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Figure  3.5  -  Block  Diagram  of  Dynamic  Wheelclimb  Computer  Simulation  Program 
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main  program  calls  a  number  of  subroutines  that  evaluate  the  first  derivatives 
of  a  set  of  simultaneous  nonlinear  state  equations,  and  then  integrates  them 
numerically.  The  program  is  written  in  FORTRAN  (IBM  extension  to  FORTRAN 
ANS) ,  and  is  compatible  with  the  IBM  G  and  H  compilers  (see  Appendix  D  for  listing) . 

The  numerical  integration  is  performed  using  DERK,  a  program  based  on  the 
sixth-order  Runge  Kutta  method,  which  is  included  in  the  IBM  library  of 
scientific  programs  [16].  DERK  finds  approximations  to  the  solution  of  a 
system  of  first-order  ordinary  differential  equations  of  the  form  x'«f (x,t) 
with  initial  conditions. 

The  derivatives  of  the  state  variables  depend  on  the  contact  geometry 
parameters  (which  are  functions  of  the  wheel set  lateral  displacement),  the 
creep  forces  (which  are  functions  of  the  wheelset  velocities,  the  contact 
geometry  and  the  normal  forces),  the  wheel/rail  contact  normal  forces,  and  the 
externally  applied  loads.  To  reduce  computation  time,  the  contact  geometry 
functions  and  creep  force/creepage  relationships  are  pre-computed  and  stored 
as  data  files  for  use  by  WHSET.  These  data  files  are  generated  by  computer 
programs  previously  developed  under  DOT  sponsorship,  known  as  WHRAILA  and 
FORCES,  and  are  widely  available  to  researchers  in  the  rail  vehicle  dynamics 
area. 

Contact  Geometry  Functions 

The  program  NHRAILA,  developed  by  R.  Heller  and  N.  Cooperrider  [17], 
solves  for  the  following  contact  geometry  parameters  as  functions  of  the 

wheelset  lateral  position:  the  wheel/rail  contact  angles  o,  the  wheelset  roll 

* 

angle  $w,  the  wheelset  c.m.  vertical  displacement  Zw,  the  rolling  radii  of 
each  wheel  ry,  the  distance  between  the  wheelset  c.m.  and  the  contact  points 
of  the  wheels  L,  and  the  lateral  radii  of  curvature  of  the  rails  R£.  WHRAILA 

*  In  this  study  the  wheelset  center  of  mass  (c.m.)  is  defined  as  the  geometric 
center  of  the  axle. 
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has  been  supplemented  by  additional  routines  which  calculate  the  values  of  the 
elliptic  integrals  a(0)  and  b(8),  the  non-dimensional  spin  parameter  x>  the 
geometric  parameter  p,  the  contact  ellipse  axes  ratio  a/b,  and  the  partial 
derivatives  defined  in  Equations  (3.4)  and  (3.5)  WHRAILA  accepts  as  inputs 
specific  wheel  and  rail  profiles,  gauge,  flange  clearance,  cant  angles  and 
rail  heights.  The  numerical  results  presented  here  are  for  new  Association  of 
American  Railroads  (AAR)  profile  wheels  running  on  rails  similar  to  132  lb  RE, 

CF  and  I  unworn  profile.  The  program  assumes  that  single  point  contact  is 
sustained  on  both  wheels  and  that  yaw  angle  effects  may  be  neglected.  Since 
these  functions  depend  on  lateral  displacement  alone,  they  do  not  have  to  be 
computed  at  each  time  step  using  WHRAILA.  Instead,  the  values  pre-computed  at 
intervals  of  0.51  mm  are  used  as  inputs  to  an  interpolation  routine.  The 
program  WHSET  uses  the  IBM  library  program  IVP,  based  on  the  Aiken’s  Lagrange 
method  of  interpolation  [16] . 

Creep  Force  Calculation 

The  second  data  file  represents  the  relationships  between  non-dimensional 
creep  forces  (f^  and  f2)  and  non-dimensional  creepages.  The  program  FORCES  [18] 
is  used  to  pre-compute  non-dimensional  creep  forces  using  the  Kalker  Simplified 
Theory.  It  takes  as  input  the  contact  ellipse  dimensions,  non-dimensional  spin 
and  creep  parameters,  and  Poisson's  ratio.  For  every  contact  ellipse  ratio  a/b 
and  non-dimensional  spin  parameter  x>  a  different  table  is  generated  for  the  creep 
coefficients  as  factions  of  the  non-dimensional  creep  parameters  in  polar  form. 

The  contact  ellipse  ratio  and  the  non-dimensional  spin  parameter  are 
determined  by  the  contact  geometry  alone,  the  latter  depending  on  the  lateral 
displacement  of  the  wheelset.  Consequently,  these  two  variables  can  be  pre¬ 
computed  for  the  full  range  of  the  lateral  motion.  A  survey  of  nunerical  results 
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over  the  displacement  range  has  been  conducted,  and  21  pairs  have  been  distinguished 
so  that  the  contact  ellipse  ratio  and  the  non-dimensional  spin  parameter  at  any 
given  displacement  are  within  10%  tolerance  of  at  least  one  of  these  21  pairs. 

The  program  FORCES  was  run  repeatedly  for  these  pairs  resulting  in  21  table  sets 
stored  as  data  files. 

At  each  lateral  displacement,  NHSET  calls  a  subroutine  to  determine  which 
of  the  21  data  files  is  to  be  used.  It  then  interpolates  the  data  in  tables  of 
this  file,  returning  the  lateral  and  longitudinal  creep  forces  as  functions  of  the 
creepages . 

Program  Outputs 

The  following  output  variables  are  generated  by  the  simulation  program: 

1.  Wheelset  lateral  displacement 

2.  Truck  lateral  displacement 

3.  Wheelset  axle  Tolling  velocity 

4 .  Wheelset  yaw  angle 

5.  Normal  force  at  each  contact  point 

6.  Derailment  quotient  for  each  wheel 

The  derailment  quotient,  the  ratio  of  lateral  to  vertical  wheel/rail 
forces,  is  the  most  frequently  used  measure  of  proneness  to  derailment  [5-10,19], 

In  the  literature  it  is  designated  by  either  (L/V)  or  (Q/P).  Following  the 
former  convention: 


l  yf2.iC05Vsinai 

'V'i  "  j  sinai-cosoi 


(3.30) 
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3.3  SIMULATION  RESULTS 

3.3.1  Definition  of  Simulation  Conditions 

The  parametric  values  for  the  simulation  are  selected  for  a  typical  high 
speed  passenger  car  described  in  [20],  and  are  sunmarized  in  Table  3.1. 

TABLE  3.1 

Simulation  Parameters 


Symbol 


Value 


Symbol 


Value 


m 

V 

1464  kg 

■t 

xx’1** 

2928  kg 

2 

680  kg-m 

I 

yy 

136  kg-m2 

u 

0.3 

e 

0.3 

E 

2xlOU  n/m2 

Fzo 

1.36xl05  n 

\ 

0 

V 

50  a/s 

2xl06  n/m 
7.4x10"*  n-m/rad 
1.17xl07  n/m 

4x10*  n-s/m 

5.84x10*  n-s/m 

2 

1.5x10  nn*$/nd 

a 

5.43x10  n-m/rad 
2.17xl02  n-m-s/rad 
0 


All  other  values  for  parameters  in  the  equations  of  motion  are  provided 
by  programs  MHRAILA  and  FORCES,  with  extensions  as  described  previously. 

Some  consents  are  needed  regarding  inclusion  of  the  truck  lateral  degree 
of  freedom.  In  many  previous  studies,  and  in  the  initial  phases  of  this 
investigation,  the  dynamics  of  the  wheelset  alone  were  simulated.  The  dynamic 
fluctuations  in  normal  force  during  derailment  simulations  that  result  from 
"wheel set -only"  models  are  unrealistically  large,  due  to  high  wheelset 


L 


accelerations.  Adding  the  truck  n ass  to  the  aodel  yields  auch  aore  natural 
behavior.  Although  the  dynamics  of  a  complete  truck  with  secondary  suspension, 
two  axles,  etc.  are  not  modeled  here,  the  present  simulations  do  exhibit 
the  important  characteristics  of  wheelsets  during  derailment. 

The  simulation  program  accepts  as  inputs  any  combination  of  forces  and 
moments  applied  to  the  truck  frame  and  wheelset,  plus  track  irregularities. 

For  clarity  in  presentation,  in  the  simulation  results  that  follow  the  only 
input  is  lateral  force  applied  to  the  truck  frame. 

3.3.2  Simulation  Results 

Representative  simulation  results  are  shown  in  Figures  3.6  and  3.7. 

Starting  with  the  wheelset  and  truck  frame  in  the  track  center  position,  the 
lateral  force  applied  to  the  truck  frame  produces  a  damped  sinusoidal  truck 
lateral  displacement.  The  wheelset  lateral  displacement  lags,  sharply  decelerating 
after  the  initial  flange  contact. 

In  the  derailment  case,  the  wheelset  first  climbs  the  rail,  falls 
slightly,  then  climbs  again  to  complete  derailment  as  the  truck  frame  dis¬ 
placement  reaches  its  second  peak.  Depending  on  initial  conditions,  the 
wheelset  may  derail  on  initial  impact  or  later,  emphasizing  the  importance  of 
the  dynamic  coupling  between  wheelset  and  truck  frame.  Even  though  the 
derailment  quotient  exceeded  the  quasisteady  limit  (2.0  for  ^w"0)  for  about 
0.05  seconds  (the  JNR  limit),  the  wheelset  does  not  derail  at  this  point. 

After  flange  contact  occurs,  the  axle  speed  drops,  as  expected  from 
quasisteady  theory.  There  is  no  fundamental  difference  between  the  derail¬ 
ment  and  non-derailment  cases. 
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Figure  3.6  -  Wheel set  aodel  simulations  for  non-derailing  (dark  line) 
and  derailing  (light  line  with  dot)  cases.  Responses  are 
to  step  input  in  lateral  force  applied  to  truck  frame. 
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The  drop  in  axle  speed  results  in  longitudinal  creep  forces  imbalance  and 
a  yaw  moment  which  steers  the  wheelset  in  the  negative  direction.  Due  to  the 
rather  stiff  primary  suspension,  the  magnitude  of  the  resulting  steering 
angles  is  quite  small  (less  than  O.S  milliradians,  or  0.03  deg).  These  small 
yaw  angles  considerably  reduce  the  coupling  between  the  longitudinal  and 
lateral  dynamics  of  the  wheelset. 

Figure  3.7  shows  the  normal  force  and  derailment  quotient  for  the  climbing 
wheel.  In  the  non-derailment  case,  the  normal  force  starts  at  a  level  close 
to  one-half  the  axle  load.  During  flange  contact,  the  normal  force  is  much 
larger  due  to  the  increasing  contact  angle.  The  rapid  deceleration  of  the 
wheelset  results  in  a  sharp  peak  in  the  normal  force.  Terms  due  to  wheelset 
lateral,  vertical,  and  roll  accelerations  all  contribute  to  the  increased 
normal  force  on  the  climbing  wheel  while  decreasing  it  on  the  opposite  wheel. 
The  initial  peak  in  normal  force  is  also  seen  in  the  L/V  response. 

In  the  derailment  case,  the  peaks  in  normal  force  and  derailment  quotient 
are  much  larger.  These  peaks  may  be  unrealistically  large,  due  to  the 
assumption  of  rigid  rails.  They  are  not  critical  to  the  simulation,  however, 
due  to  their  extremely  short  duration.  Peaks  shown  after  initial  flange  contact 
are  due  to  large  values  of  partial  derivative  terms  in  Equations  (3.10)  and 
(3.11). 

Just  prior  to  derailment,  the  normal  force  goes  to  zero  for  a  short  time. 
Physically,  the  wheel  is  beginning  to  lose  contact  while  the  contact  angles 
decreases.  The  large  fluctuations  in  normal  force  and  loss  of  contact  prior  to 
derailment  have  been  confirmed  experimentally  in  tests  at  JNR  [10]. 
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3.3.3  Discussion  of  Modeling  Effects 

The  computer  program  was  altered  to  study  the  effect  of  model  complexity 
on  the  simulation  results.  The  first  simplification  was  to  set  the  axle  speed 
to  V/rQ.  Exclusion  of  this  degree  of  freedom  was  found  to  have  almost  no 
effect  on  the  simulation  of  the  other  dynamic  variables  as  a  result  of  the 
stiff  suspension  and  the  negligible  steering  effects.  As  a  consequence,  the 
value  of  the  truck  yaw  angle  iji  is  very  important  to  derailment  safety,  since 
the  derailment  limit  falls  from  2.0  to  0.8  as  increases  from  0  to  SO  milli- 
radians.  This  result,  predicted  from  quasisteady  derailment  theory,  is  verified 
by  the  dynamic  simulation  model.  Passenger  trucks  with  softer  yaw  suspensions 
and  freight  trucks  may  exhibit  different  behavior. 

The  acceleration  terms  in  Equations  (3.11)  through  (3.26)  are  significant 
in  terms  of  their  effects  on  wheelset  motion,  derailment  quotient,  and  derail¬ 
ment  proneness.  In  contrast,  gyroscopic  terms  are  not  significant  and  may  be 
omitted. 

The  simulation  model  presented  in  this  report  represents  a  successful 
extension  of  contemporary  understanding  on  wheel/rail  contact  phenomena  and 
wheelset  dynamics  to  the  wheelclimb  derailment  regime.  Through  the  use  of 
this  detailed  model,  it  has  been  determined  that  the  following  elements  are 
necessary  for  simulation  of  wheelclimb  derailment  of  a  typical  high  speed 
passenger  car  with  a  stiff  yaw  suspension: 

1.  Wheelset  and  truck  lateral  displacement  degrees  of  freedom. 

2.  Wheel/rail  contact  geometry  parameters  determined  as  functions  of 
wheel  and  rail  profiles. 

3.  Creep  forces  determined  by  Kalker  Simplified  Theory. 


L 


3-25 


4.  Wheel/rail  normal  forces  computed  from  kinematic  constraints  resulting 
from  continuous  two-wheel,  single-point  contact,  including  acceleration 
effects. 

5.  Creepages  determined  by  wheelset  kinematics. 

6.  Truck  yaw  angle. 

For  freight  trucks  or  passenger  trucks  with  soft  primary  suspensions, 
wheelset  yaw  angle  and  axle  rotational  speed  should  be  included  as  additional 
degrees  of  freedom.  Gyroscopic  terms  are  not  critical  to  wheelclimb  simulation. 

The  reduction  in  necessary  degrees  of  freedom  associated  with  stiff  yaw 
suspension  trucks  results  in  significant  reductions  in  computation  time, 
which  will  be  particularly  important  in  multiple  axle,  complete  vehicle 
simulations. 
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3.4  Experimental  Validation  and  Usefulness  of  Model 
3.4.1  Results  of  Phase  One  Experiments 

To  validate  the  theoretical  analysis  and  simulation  model,  a  series  of 
experiments  on  dynamic  wheelclimb  was  conducted  during  Phase  One  of  this 
research  program.  The  apparatus  used  is  described  in  Appendix  B,  being  identical 
to  that  used  for  study  of  quasisteady  wheelclimb.  Forces  acting  on  the  wheel- 
set  were  measured  with  the  six-component  strain  gauge  balance,  so  that  only  axle 
L/V  ratios  could  be  measured.  The  only  addition  to  the  apparatus  used  in  the 
quasisteady  wheelclimb  tests  was  a  photocell  actuated,  pneumatic  servo  to  apply 
step  or  pulse  inputs  in  lateral  load  to  the  wheelset  at  designated  locations  along 
the  track. 

To  simulate  the  mechanics  of  the  wheelset  test  apparatus,  it  was  necessary 
to  modify  the  equations  of  motion  derived  in  Section  3.1.  These  changes  were 
required  to  properly  represent  the  kinematics  and  dynamics  of  the  rigid  body 
elements  that  comprise  the  gimbel  and  linkage  system  described  in  Appendix  B. 

Since  the  modified  equations  apply  only  to  this  unique  apparatus  and  are  not  of 
general  interest,  they  are  not  included  in  this  report  (for  details,  refer  to 
[15]). 

Representative  simulated  and  measured  responses  for  non-derailing  and 
derailing  cases  are  shown  in  Figures  3.8  and  3.9.  The  simulations  include  the 
dynamics  of  the  pneumatic  actuator  used  to  apply  the  lateral  force,  which  acts 
as  a  first-order  lag  with  a  time  constant  of  125  ms.  The  agreement  between  pre¬ 
dicted  and  measured  lateral  displacements  for  non-derailing  and  derailing  cases 
is  good,  although  for  the  latter  case  the  predicted  derailment  occurs  more 
rapidly  than  that  measured. 
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Significant  differences  between  predicted  and  measured  time  histories 
of  the  L/V  ratio  are  evident  in  the  figures  shown.  Specifically,  the  pre¬ 
dicted  responses  characteristically  exhibit  higher  maximum  values  and  shorter 
time  durations  than  those  measured.  To  resolve  this  discrepancy,  further 
research  conducted  under  Phase  Two  of  this  study  was  directed  towards  refining 
the  analytical  model  and  improving  the  force  measuring  system  used  in  the 
experiments.  The  results  of  Phase  Two  experiments  are  presented  in  the 
following  section. 

3.4.2  Results  of  Phase  Two  Experiments 

At  the  onset  of  Phase  Two,  there  were  two  principal  factors  considered 
as  having  the  potential  for  influencing  the  measured  wheel  force  signals: 

1.  The  strain  gauge  balance  signals  measured  in  Phase  One  did  not 

resolve,  in  an  unambiguous  manner,  the  various  external  and  internal 
forces  acting  during  the  dynamic  wheelclimb  experiment.  Under  quasi¬ 
steady  conditions,  the  external  forces,  in  both  vertical  and  lateral 
directions,  approximately  balance  wheel/rail  forces  acting  at  the 
contact  points.  Since  the  strain  gauge  balance  was  mounted  in  series 
between  the  external  and  wheel/rail  forces,  the  strain  gauge  balance 
was  capable  of  sensing  the  desired  variables  under  quasisteady  conditions. 
Further,  under  quasisteady  conditions,  a  simple  functional  relation 
exists  relating  axle  and  wheel  L/V  ratios,  so  that  resolution  of 
independent  wheel  ratios  is  not  necessary.  Under  dynamic  wheelclimb 
conditions,  none  of  the  above  apply,  so  that  it  is  desirable  to  have 
independent  sensing  of  individual  wheel  loads,  external  forces,  and 
inertial  forces  (i.e.  accelerations). 
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2.  The  bandwidth  of  the  instrumentation  system  used  in  Phase  One  was  25  Hz. 
The  25  Hz  cutoff  frequency  was  determined  to  be  necessary  to  prevent 
aliasing  of  signals  passed  through  the  multiplexing  circuits  prior 
to  recording  on  the  FM  tape  recorder.  The  L/V  ratio  pulses  predicted 
from  the  digital  simulation  contain  significant  high  frequency  content 
that  would  be  attenuated  by  a  25  Hz  lowpass  filter. 

To  eliminate  these  limitations  in  the  dynamic  wheelclimb  experiments,  a 
new  set  of  apparatus  for  force  measurement  was  designed,  fabricated,  and  carefully 
calibrated.  Digital  simulations  were  performed  which  passed  the  predicted  wheel 
lateral  and  vertical  force  signals  through  lowpass  filters,  to  determine  the 
bandwidth  necessary  for  recording  of  these  signals  in  an  unaltered  state.  The 
simulations  clearly  showed  that  the  25  Hz  filter  used  in  Phase  One  significantly 
attenuated  and  widened  the  L/V  ratio  pulse  to  qualitatively  match  that  measured 
in  Phase  One.  These  simulations  identified  filter  bandwidth  as  the  principal 
cause  of  the  Phase  One  discrepancy  between  theory  and  test.  Further  simulations 
showed  that  a  cutoff  frequency  of  100  Hz  was  sufficient  for  recording  wheel  force 
signals  without  attenuation. 

An  instrumented  wheelset  was  designed  to  measure  lateral  and  vertical  forces 
at  each  wheel/rail  contact  point.  The  ASEA/SJ  configurations  for  lateral  and 
vertical  bridges  were  selected.  Strain  gauge  locations  were  determined  after  map¬ 
ping  of  the  strain  fields  on  the  wheel  plate  surfaces,  to  maximize  signal  output 
with  minimum  ripple,  load  point  sensitivity,  and  cross  talk.  External  forces  were 
measured  separately  with  load  cells,  with  the  strain  gauge  balance  used  in  Phase  One 
installed  to  provide  capability  for  Phase  One/Phase  Two  data  comparisons.  Details 
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of  the  apparatus  design  and  calibration  are  given  in  Appendix  C. 

Representative  measurements  of  signals  during  derailment  using  the 
improved  instrumentation  system  are  presented  in  Figures  3.10  through  3.17. 

At  the  beginning  of  the  test,  the  wheelset  is  positioned  in  tread  contact 
near  the  opposite  (left)  rail  using  a  special  fixture.  At  a  designated 
track  location,  a  photocell  senses  a  light,  activating  the  lateral  force 
actuator.  The  lateral  force  rises  exponentially,  as  shown  in  Figure  3.10. 

The  wheelset  moved  laterally  in  response  to  this  input  until  it  makes  flange 
contact  on  the  right  rail,  where  it  pauses  for  about  22  msec.  After  the 
wheel  has  climbed  the  rail,  lateral  displacement  increases  rapidly  to  complete 
>  derailment,  as  shown  in  Figure  3.11.  As  the  wheelset  approaches  one  rail, 

a  strong  yaw  moment  is  produced  by  the  longitudinal  creep  forces,  causing  a 
steering  action  in  the  negative  direction,  as  shown  in  Figure  3.12.  As 
demonstrated  in  both  figures,  the  agreement  with  simulation  is  excellent. 

Figures  3.13  through  3.15  show  the  forces  acting  at  the  wheel/rail  interfaces. 
A  large  pulse  in  lateral  force  on  the  derailing  wheel  is  produced  during  the 
period  that  the  wheel  is  in  flange  contact  (about  40  msec).  This  short  duration 
is  consistent  with  theory,  and  clearly  shows  that  the  long  duration  (order  of  100 
msec.)  measured  in  Phase  One  resulted  from  the  25  Ht  filters.  The  lateral  forces 
measured  on  the  non-derailing  wheel  result  from  creepage,  and  are  consistent  with 
a  friction  coefficient  of  0.3.  The  vertical  force  on  the  derailing  wheel 
(Figure  3.15)  shows  some  fluctations  from  its  nominal  value  of  22  n.  The 
excursions  at  the  flange  contact  point  are  consistent  with  model  predictions, 
resulting  from  accelerations  in  the  vertical  direction  caused  by  contact  geometry 
constraints.  Other  fluctuations  in  vertical  force  result  from  track  disturbances 
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in  the  vertical  direction  that  were  not  included  in  the  simulation  model  nor 
measured  in  the  experiment. 

The  lateral  acceleration  of  the  wheelset  is  shown  in  Figure  3.16.  The 
maximum  deceleration  is  about  -0.2  g,  which  corresponds  to  an  inertial  force  of 
about  -22  n.  When  compared  to  the  lateral  force  measured  at  the  derailing  wheel 
(about  75  n.  at  the  same  time),  the  inertial  force  is  significant.  This  indicates 
that  measurement  of  wheel/rail  loads  at  locations  as  close  as  possible  to  the 
contact  points  (i.e.  using  instrumented  wheelsets)  is  superior  to  measurement 
at  an  interior  point  (i.e.  using  strain  gage  balance). 

Finally,  the  predicted  and  measured  L/V  ratios  are  shown  in  Figure  3.17. 

The  amplitude  and  duration  of  the  primary  and  secondary  peaks  in  the  response 
are  predicted  well  by  the  simulation.  The  first  peak  results  from  the  large  am¬ 
plitude  pulse  in  lateral  force  during  flange  contact;  the  second  peak  results 
from  the  dip  in  vertical  force  resulting  from  wheelset  vertical  acceleration. 

The  experiments  performed  in  Phase  Two  demonstrate  that  the  dynamic  wheel - 
climb  model  presented  in  this  report  is  capable  of  predicting  both  wheelset 
response  to  external  force  inputs  and  the  resulting  wheel/rail  contact  forces. 

This  validated  model  should  prove  useful  in  the  future  in  the  evaluation  of 
vehicle  derailment  safety  and  synthesis  of  new  derailment  criteria. 
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Figure  4.16  Generic  Lateral  Velocity  Responses  for  Process  B  Dynamic 

Wheelclimb  Tests.  Double  Peak  Indicates  Marginal  Derailment 
Condition. 
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Figure  3.11  -  Measured  and  simulated  wheel set  lateral  displacement. 
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Figure  3.12  >  Measured  end  simulated  wheelset  yaw  angle 
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4.  EVALUATION  OF  DERAILMENT  CRITERIA 
4-1  Commonly  Used  Criteria 

Although  a  larger  set  of  variables  are  available  from  the  experiments 
and  simulations  for  use  in  synthesizing  criteria  for  wheelclimb  derailment, 
in  this  report  attention  is  restricted  to  those  criteria  which  make  use  of 
the  time  histories  of  wheel  loads  only.  The  objective  here  is  to  evaluate 
the  hypothesis  that  sufficient  information  exists  in  the  wheel  load  measure¬ 
ments  alone  to  predict  derailment  safety  reliably. 

The  basic  information  available  in  these  pulselike  signals  may  be 
distilled  into  the  following  quantities:  peak  value,  integrated  value, 
average  value  over  pulse  duration,  pulse  shape,  duration  above  threshold, 
etc.  Table  4.1  is  a  summary  of  candidate  criteria  which  are  of  the  general 
form  of  a  load  amplitude  versus  time  duration.  To  test  the  hypothesis  a 
wide  selection  of  candidate  criteria  was  explored.  The  widely  used  JNR 
Criterion  is  of  Type  4  in  the  table;  duration  t^  is  defined  to  be  1.5  times 
the  period  during  which  the  pulse  exceeds  50%  of  its  peak  value  [30], 

Criteria  1  through  5  yield  a  single  data  point  for  each  event; 

Criteria  6  yields  a  locus  of  points  as  the  threshold  value  is  swept  from 
zero  to  the  pulse  peak.  All  criteria  may  be  applied  to  either  L/V  ratios 
or  lateral  forces  alone.  Criteria  involving  peak  values  may  lack  physical 
significance  and  are  vulnerable  to  erroneous  conclusions  drawn  from  noisy 
instrumentation  or  wheel  loads  resulting  from  high  frequency  vibration. 
Integrated  lateral  force  has  units  of  change  in  momentum,  although  it  does 
not  represent  the  time  change  in  momentum  since  only  the  wheel/rail  reaction 


TABLE  4.1 


in  a 


forces  are  being  included.  By  dividing  the  integrated  value  by  the  time 
duration,  the  magnitude  averaged  over  the  pulse  is  obtained.  Criteria 
4  and  5  may  be  used  to  define  pulse  duration  when  the  mean  value  of  the 
signal  before  and  after  the  pulse  is  above  zero. 

The  sensitivity  of  each  criterion  to  changes  in  pulse  area,  shape, 
and  duration  is  shown  in  Table  4.1.  in  this  context  sensitivity  refers 
to  the  ability  of  a  criterion  to  discriminate  between  pulses  on  the  basis 
of  area,  duration,  or  shape  by  yielding  a  different  (x,y)  value  or  locus 
in  the  graph.  For  example,  if  a  criterion  sensitive  to  area  were  found 
to  be  a  good  derailment  predictor,  then  derailment  could  be  correlated 
with  momentum.  Many  additional  cr  ? • ia  can  be  proposed,  but  they 
generally  will  perform  in  a  manner  very  similar  to  those  in  Table  4.1. 

4.2  Experimental  Evaluation  of  the  JN’R  Criterion 

The  recorded  experimental  wheel  load  and  lateral  displacement  time 
histories  were  digitized  and  analyzed  using  computer  subroutines  which 
implemented  the  various  criteria. 

For  Process  B,  112  derailment  events  (occurrence  of  large  L/V  ratios 
or  lateral  displacements)  were  examined,  distributed  in  the  test  matrix 
in  Table  4.2.  Cases  were  selected  deliberately  to  be  near  the  derailment 
limit,  and  do  not  reflect  the  distribution  of  their  occurrence  in  the 
field.  A  special  feature  of  these  experiments  is  the  large  nuftber  of 
marginal  and  complete  derailments,  which  would  not  generally  be  available 


from  field  tests. 
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TABLE  4.2 

TEST  MATRIX  FOR  DYNAMIC  WHEELCLIMB  EXPERIMENTS 


TEST  MATRIX 


Yaw  Angle 


3° 

0° 

-3° 

Total 

No  Derailment 

id 

19 

1 

30 

Marginal  Derailment 

9 

13 

17 

39 

Complete  Derailment 

8 

21 

14 

43 
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Figure  4.1  illustrates  the  application  of  JNR  criterion  to  a  typical 
L/V  ratio  pulse  using  the  algorithm  shown  in  Figure  4.2 »  with  the  results 
of  similar  application  to  the  events  in  the  test  matrix  shown  in  Figure  4.3. 
The  JNR  Criterion  is  modified  using  Eqs.  (2.4)  and  (2.5)  to  convert  wheel 
L/V  ratio  to  axle  L/V  ratio.  It  is  clear  from  these  results  that  not  only 
does  the  JNR  criterion  not  predict  derailment,  but  that  no  criterion  based 
solely  on  L/V  ratio  and  time  duration  could  separate  the  derailment,  marginal 
derailment,  and  nonderailment  regions.  The  data  in  Figure  4.3  includes 
measurements  at  all  yaw  angles;  from  the  quasisteady  derailment  results 
strong  influence  of  yaw  angle  could  be  expected  by  rescaling  the  data  and 
JNR  Criterion  by  the  quasisteady  limits  shown  in  Figure  4.2  (0.4  at  *3  deg., 

0.6  at  0  deg.,  3.1  at  -3  deg.),  a  normalized  version  of  the  JNR  Criterion 
results,  given  in  Figure  4.4.  The  JNR  Criterion  is  still  an  inadequate 
derailment  predictor  using  this  modification. 

4 . 3  Experimental  Evaluation  of  Alternative  Wheel .Load  Criterion 

The  alternative  wheelload  derailment  criteria  are  evaluated  using 
the  experimental  data  base  by  computing  the  following  quantities: 

TAT  *  time  above  threshold 

PEAK  «  peak  value  of  L/V  ratio 

AREA  -  LA  -  /  (rr)  dt 
TAT  v 

MOMENTUM  «  /TAT  L  dt 

(AREA- L/V) /TAT  -  ^  /TAT  (£)  dt 

NORMALIZED  MOMENTUM  ■  yjy  /?AT  L  dt 

TJNR  «  tj  in  Figure  4.1 


Figure  4.1  -  Application  of  JNR  Derailment  Criterion  to  Typical  L/V  Pulse 
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Figure  4.2  •  Algorithm  for  Derailment  Criteria  Evaluation 
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4.3.b.  Compariaon  of  Dynamic  Wheel  climb  Experimental  Data  with 
JNR  Criterion  for  Marginal  Derailment  Caie*. 
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Figure  4.3.C.  Comparison  of  Dynamic  Wheelclimb  Experimental  Data  with 
JNR  Criterion  for  -Complete  Derailment  Cases. 
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Figure  4.4.b.  Comparison  of  Dynamic  Wheelclimb  Experimental  Data  with  JNR 
Criterion  with  L/V  Ratios  Rescaled  by  the  Quasisteady  - 
Derailment  Limit  at  Each  Yaw  Angle  for  Marginal  Derailment  Cases 
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Figure  4.4.c. 
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Comparison  of  Dynamic  Wheelclimb  Experimental  Data  with  JNR 
Criterion  with  L/V  Ratios  Rescaled  by  the  Quasisteady 
Derailment  Limit  at  Each  Yaw  Angle  for  Complete  Derailment  Cases 
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The  results  of  these  computations  are  tabulated  in  Appendix  A.  In 
the  following  paragraphs  the  results  of  application  of  the  various  criteria 
•re  plotted  and  discussed. 

Type  l  -  Peak  value  of  L/V  versus  TAT 

Figure  4.5  shows  data  for  the  three  derailment  conditions  plotted 
in  absolute  and  scaled  formats  (using  L/V  ratio  scale  factors  for  each  yaw 
angle).  As  is  the  case  using  the  JNR  criterion,  the  spread  of  peak  values 
for  each  derailment  case  covers  almost  a  factor  of  ten  in  the  absolute 
format  and  a  factor  of  six  in  the  scaled  format.  No  significant  time 
duration  dependence  is  evident. 

Type  2  -  Integrated  Values  versus  TAT 

Figure  4.6  and  4.7  show  data  for  L/V  ratio  and  lateral  force  integrated 
over  the  pulse  duration,  the  latter  yielding  a  quantity  with  units  of 
momentum.  However,  as  previously  discussed,  the  momentum  quantity  plotted 
here  is  not  the  momentum  of  the  wheelset,  but  rather  the  change  in 
momentum  due  to  action  of  the  contact  forces  alone.  In  both  figures  all 
data  follow  a  distinct  trend  of  increasing  integrated  values  with  increasing 
pulse  durations.  This  result  is  consistent  with  the  observation  that  the 
longer  pulses  are  rather  rectangular  in  shape.  Neither  criteria  discriminates 
successfully  among  the  three  derailment  conditions. 

Type  3  -  normalized  Integrated  Values  versus  TAT 

The  proportionality  between  integrated  values  in  Type  2  criteria  and 
time  duration  is  compensated  by  dividing  the  integrated  value  by  the  time 
duration.  The  resulting  data  are  then  the  average  values  computed  over  the 


Figure  4.6  Absolute  and  Scaled  Derailment  Test  Data:  Integrated  L/V  Ratio  Versus  Time 
Above  Threshold. 
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duration  of  the  pulse.  Figures  4.8  and  4.9  show  NORMALIZED  AREA-L/V  and 
MOMENTUM,  respectively.  The  results  are  again  very  similar  to  those 
for  the  JNR  and  Type  1  criteria. 

Type  4  -  Peak  Value  versus  TNP 

The  JNR  Criterion  is  of  this  type,  with  results  presented  previously. 
Type  S  -  Integrated  Value  versus  TNP 

Same  results  as  for  Type  2,  with  data  shifted  to  the  left  in  each 
plot,  since  TNP  is  always  less  than  TAT.  Since  no  significant  influence 
of  time  duration  has  been  found  in  the  test  results.  Types  2  and  3 
criteria  are  equivalent. 

Type  6  -  Threshold  Value  versus  Exceedance  Tine 

This  type  of  criterion  is  applied  by  sliding  the  threshold  value, 
which  is  plotted  as  the  ordinate,  from  zero  to  the  peak  value  of  the 
signal  and  computing  the  time  that  the  pulse  exceeds  the  threshold.  In 
this  manner,  the  pulse  is  mapped  into  the  threshold  -  exceedance  time 
plane  as  a  continuous  curve  rather  than  a  single  point.  Since  the  continuous 
curve  retains  much  of  the  information  content  of  the  pulse,  in  contrast 
to  criteria  which  reduce  the  pulse  information  to  a  single  point,  this 
procedure  may  have  more  potential  for  derailment  condition  discrimination. 

The  exceedance  curves  plotted  are  for  single  pulses,  and  should  not  be 
confused  with  statistical  plots  for  continuous  running  time  histories 
including  many  pulses  and  intervening  periods. 


Figure  4.10  Threshold  Level  Versus  Exceedance  Time  Plots  for  Derailment  Test  Pulses. 
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The  results  of  this  procedure  applied  to  L/V  ratio  pulses  are  shown 
in  Figure  4.10.  The  shapes  of  the  curves  plotted  are  generally  similar 
to  one  another,  indicative  of  pulse  shapes  ranging  from  half  sine  waves 

1 

to  rectangles.  The  curves  for  the  three  derailment  conditions  generally 
overlap  in  threshold  values  and  exceedance  times,  with  no  apparent  divisions 
among  them. 

4.4  Results  from  Variable  Yaw  Angle  Tests 

Tests  were  conducted  with  the  yaw  degree  of  freedom  unlocked  and 
a  soft  linear  spring  yaw  suspension  installed.  These  experiments  for 
Process  C  permitted  self  steering  by  the  wheelset  under  flange  contact 
conditions  and  hunting  phenomena  at  high  speed.  The  data  recorded  were 
most  interesting,  but  in  retrospect  the  number  of  data  points  accumulated 
was  too  small  to  yield  statistically  meaningful  results.  The  section  that 
follows  discusses  the  implications  of  the  recorded  observations  as  examples 
of  the  complexity  of  the  dynamics  of  Process  C,  but  at  this  time  conclusions 
as  definite  as  those  for  Processes  A  and  B  cannot  be  drawn. 

The  four  generic  types  of  responses  recorded  for  the  wheelset  are 
shown  in  Figures  4.11  through  4,14.  In  Figure  4.11  the  normal  response 
of  the  wheelset  at  sub-critical  velocities  to  track  lateral  random  mis¬ 
alignments  is  shown.  In  the  absence  of  flange  contact  axle  L/V  ratios 
are  small  (less  than  0.05  ),  lateral  displacements  relative  to  the^ track 
less  than  ±0.25  mm,  and  yaw  angles  less  than  10.S0  mm.  In  Figure  4.12  limit 
cycle  hunting  of  the  wheelset  is  shown  for  velocities  in  excess  of  6  m/s 
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nod  el  speed.  The  approximate  ranges  are  ±0.4  for  axle  L/V  ratios.  —  0.5  min 
for  lateral  displacements  relative  to  the  track,  and  -  1*5°  for  yaw  angles. 
Derailments  did  not  occur  under  these  conditions. 

In  Figures  4.13  and  4.14  the  responses  to  step  inputs  in  lateral  force 
are  shown.  The  tests  were  conducted  in  the  same  manner  as  for  Process  B 
except  that  the  rubber  retaining  wheels  preventing  derailment  on  the  left 
rail  were  removed.  Axle  L/V  ratios  up  to  1.2  were  recorded,  with  yaw  angles 
in  the  range  of  *2°  to  -6°,  and  lateral  displacements  including 

full  derailment.  This  type  of  response  occurred  up  to  3  m/s  model  speed. 
Above  3  m/s  but  below  the  apparent  critical  speed  of  6  m/s  a  limit  cycle 
hunting  motion  was  initiated  after  release  of  the  applied  lateral  force. 

This  response  is  shown  in  Figure  4.14.  The  wheelset  may  or  may.  not  derail 
due  to  the  initial  application  of  the  lateral  force,  but  once  the  force 
was  removed  and  hunting  started  no  further  derailments  were  observed. 
Detailed  examination  of  this  limit  cycle  showed  that  it  is  identical  to 
the  one  found  above  6  m/s  in  frequency  and  trajectory  in  the  phase  plane. 

The  stability  of  the  limit  cycle  behavior  is  shown  in  Figure  4.15,  with 
a  stable  limit  cycle  existing  for  speeds  below  the  critical  velocity, 
contrary  to  limit  cycle  behavior  predicted  in  several  recent  papers.  The 
production  of  a  limit  cycle  after  a  large  lateral  force  is  removed  is  of 
practical  significance,  since  this  could  simulate  conditions  at  the  exit 
from  a  curve.  Additional  discussion  of  this  behavior,  known  as  a  Hopf 
bifurcation  is  given  in  [2,4]. 

The  most  interesting  results  from  a  derailment  prediction  point  of 
view  are  those  typified  by  Figure  4.13.  For  the  same  applied  lateral 
force  both  derailments  and  nonderailments  occur. 
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Figure  4.13  Two  Hegree  of  Freedom  Wheelset  Response  to  Step  Lateral  Force  Inputs  in  Addition 
to  Random  Track  Alignment. 


Stability  of  Limit  Cycle  behavior  as  Function  of  Forward  Speed  and 


I 


4-29 


4  Close  examination  of  the  responses  shows  several  key  phenomena: 

1)  The  yaw  moment  acting  on  the  wheel  set  when  in  flange  contact 
causes  the  axle  to  be  steered  towards  the  track  center;  i.e. 

I  a  positive  lateral  displacement  leads  to  a  negative  yaw  rate. 

The  yaw  moment  acts  to  inhibit  the  derailment  process,  since  as 
the  yaw  angle  becomes  more  negative,  the  lateral  creep  forces 
increase  in  the  negative  direction. 

2)  The  key  determinant  of  derailment  is  the  yaw  angle  at  the  time 

of  initial  flange  contact.  In  Figure  4.13,  the  yaw  angles  at  contact 
for  the  derailing  cases  are  about  -1  degrees,  while  for  the 
nonderailing  case  it  is  about  -3  degrees.  The  value  of  'he 
yaw  angle  at  flange  contact  is  a  random  variable,  determined 
by  the  wheelset  motion  prior  to  application  of  the  lateral  force. 
While  insufficient  data  was  obtained  to  evaluate  derailment  criteria 
quantitatively,  the  observations  described  above  provide  a  qualitative 
indication  that  criteria  which  do  not  explictly  or  implictly  account  for 
yaw  angle  effects  will  not  be  successful  predictors  of  derailment  safety. 

4.5  Discussion  of  Experimental  Results  for  Dynamic  Wheclclimb  Processes 

The  experimental  data  presented  in  this  section  shows  that  a  variety 
of  time-duration  dependent  wheel  load  derailment  criteria  are  unsuccessful 
in  predicting  derailment  safety.  The  nonderailment,  marginal  derailment, 
and  complete  derailment  cases  are  represented  by  these  criteria  with  data 
points  in  overlapping  regions  spread  generally  over  an  order  of  magnitude. 
Defining  derailment  safety  criteria  below  the  minimum  levels  at  which 
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derailment  occurred  would  be  overly  conservative,  possibly  excluding 
vehicles,  operating  conditions,  and  track  maintenance  standards  that 
could  be  demonstrated  to  be  safe  with  more  discriminating  criteria. 

Even  if  these  minimum  levels  were  to  be  used,  they  do  not  establish 
a  lower  bound  for  derailment,  since  lower  data  points  possibly  could 
be  measured. 

The  data  presented  do  not  indicate  that  a  time-duration  dependence 
should  be  included  in  derailment  criteria.  The  time  durations  measured 

are  in  the  range  of  20  ms  to  50  ms,  somewhat  above  the 
range  of  durations  in  JNR  data  [7].  It  is  possible  that  under  different 
experimental  conditions  impulsive  wheel  loads  with  very  short  (less  than 
10  ms)  time  durations  and  high  amplitudes  could  be  measured  that 
would  justify  a  time-duration  dependence  in  wheel  load. 

4.6  Derailment  Diagnostics 

In  the  scale  model  test  program  reported  here  it  was  equally  feasible 
to  measure  responses  that  did  or  did  not  involve  derailment.  In  full 
scale  tests  of  vehicle  safety,  in  general  it  will  not  be  feasible  to  derail 
vehicles.  This  presents  a  dilemna  analagous  to  nondestructive  material 
testing,  being  that  if  the  test  conditions  are  not  severe  enough,  no  data 
relevant  to  safety  limits  may  be  obtained.  It  is  most  desirable  to  be  able 
to  detect  the  approach  to  a  derailment  limit  without  exceeding  it. 

The  measured  responses  shown  in  Figure  4.16  demonstrate  such  a  diagnostic 
technique.  By  high-pass  filtering  and  integrating  the  signal  from  an 


accelerometer  attached  to  the  wheelset,  the  lateral  velocity  of  the  wheel  - 
set  is  obtained.  The  single  peak  shown  in  the  left  half  of  Figure  4.16 
is  generic  to  nonderailments,  while  the  double  peak  on  the  right  half  is 
generic  to  marginal  derailments.  The  physical  explanations  for  these 
responses  are  straightforward.  At  initial  flange  contact,  the  rapid 
increase  in  contact  angle  rapidly  decelerates  the  wheelset  in  the  lateral 
direction,  resulting  in  the  first  peak.  In  a  nonderailment  the  velocity 
returns  to  zero,  the  derailment  process  stops  with  only  a  single  peak. 

In  marginal  and  complete  derailments  the  velocity  does  not  return  to 
zero  after  the  first  peak.  In  complete  derailments  after  the  deceleration 
associated  with  the  first  peak  occurs,  the  wheelset  again  accelerates. 

In  this  case  the  velocity  increases  until  the  wheelset  leaves  the  track. 

In  marginal  derailments  the  lateral  forces  applied  to  the  wheelset  reverse 
direction  before  the  complete  derailment,  centering  the  wheelset  on  the 
track.  This  last  second  rescue  causes  the  second  peak. 

This  phenomenon  may  be  readily  exploited  to  achieve  a  simple  diagnostic 
indicator  of  impending  derailment  conditions.  Since  the  number  of  peaks, 
rather  than  absolute  values  of  signals,  are  involved,  simple  logic  circuits 
or  microprocessor  computer  programs  could  test  signals  in  real  time.  For 
example,  the  forward  velocity  of  a  locomotive  under  test  could  be  increased 
on  succeeding  runs  over  a  perturbed  track  section  until  double  peaks  are 
detected.  At  this  point  this  test  sequence  could  be  terminated,  with  the 
knowledge  that  a  derailment  safety  limit  had  been  found. 


Figure  4.16  Generic  Lateral  Velocity  Responses  for  Process  B  Dynamic 

Wheelclimb  Tests.  Double  Peak  Indicates  Marginal  Derailment 
Condition. 
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While  more  research  is  required  to  establish  and  validate  diagnostic 
tools  applicable  to  full  scale  vehicles,  the  initial  success  of  this  technique 
found  from  single  wheelset  experiments  is  most  encouraging. 
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5.  CONCLUSIONS 

5.1  Conclusions 

The  experimental  and  analytical  studies  of  wheelclimb  derailment 
presented  in  this  report  result  in  the  following  conclusions  regarding 
derailment  criteria: 

1)  Wheelclimb  derailment  criteria  based  on  quasisteady  theory  are 
adequate  for  derailment  prediction  under  conditons  of  negligible 
lateral  incident  velocity  and  constant  positive  yaw  angle.  At 
negative  yaw  angles  derailment  occurs  at  L/V  ratios  somewhat 
below  the  predicted  limits. 

2)  Application  of  the  JNR  and  other  time-duration  dependent  derailment 
criterion  for  nonderailment,  marginal  derailment,  and  complete 
derailment  cases  are  each  spread  over  an  order  of  magnitude  in 

L/V  ratio  amplitude  and  time  duration.  None  of  the  criterion  tested 
could  distinguish  between  safe  and  unsafe  conditions.  Defining 
a  safety  criterion  below  the  derailment  data  would  be  overly 
conservative,  possibly  excluding  vehicles,  operating  conditions, 
and  track  and  maintenance  standards  that  could  be  demonstrated  to 
be  safe  with  a  more  discriminating  criterion.  Furthermore,  the  data 
provided  do  not  necessarily  establish  a  lower  bound  for  derailment, 
since  lower  points  possibly  could  be  measured. 

3)  The  analytical  models  for  dynamic  wheelclimb  yield  accurate 

predictions  of  wheelset  response  to  external  force  inputs, 

in  terms  of  wheelset  motions  during  derailment  and  wheel/rail 

« 

interaction  forces  (i.e.  L/V  ratios). 
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4)  Evidence  has  been  found  that  derailment  criteria  employing  variables 
measured  in  addition  to  wheel  loads  may  be  successful  in  predicting 
derailment  safety,  and  that  diagnostic  criteria  may  be  developed 
for  warning  of  impending  derailment. 

5.2  Recommendation  for  Future  Research 

To  achieve  the  objective  of  defining  a  reliable  measure  of  derailment  safety, 
the  following  would  be  useful  as  topics  of  future  research: 

a)  New  wheel climb  derailment  criteria  should  be  developed  and  validated 
that  include  variables  in  addition  to  wheel  loads,  such  as  lateral 
velocity  and  yaw  angle,  that  are  readily  measured  under  full  scale 
test  conditions.  Such  multivariable  criteria  should  be  better 
indicators  of  derailment  safety. 

b)  Criteria  for  the  wheelclimb  process  should  be  combined  with 
criteria  for  other  derailment  modes,  such  as  gauge  spreading 

and  rail  rollover,  to  yield  a  comprehensive  safety  measure.  Track 
stiffness  parameters  would  be  key  variables  in  the  comprehensive 
criteria  formulated. 

c)  The  results  of  analytical  and  experimental  studies  of  derailment 
of  single  wheelsets  should  be  extended  to  complete  trucks  and 
vehicles.  It  is  very  important  to  establish  the  degree  to  which 
single  wheelset  criteria  may  be  applied  directly  to  complete 
vehicle  configurations. 
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In  addition  to  the  above,  the  concept  developed  in  this  study  of 
derailment  diagnostics  should  be  pursued.  This  technique  provides  a  means 
for  detecting  marginally  safe  conditions  in  full  scale  tests,  so  that 
safety-related  phenomena  may  be  measured  without  actual  derailments  being 
required  in  the  test  plan.  Scale  model  experiments  will  continue  to  be 
useful  to  the  study  of  fundamental  derailment  processes,  due  to  the  greater 
control  of  test  conditions  possible  and  the  relative  ease  of  study  of  the  full 
range  of  derailment  conditions. 
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Appendix  A 


SUMMARY  OF  DERAILMENT  TEST 
DATA  AND  COMPUTED  INDICES 
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APPENDIX  B 


DESCRIPTION  OF  PHASE  ONE  APPARATUS  AND  WHEELCLIMB  EXPERIMENTS 

B.l  Description  of  Wheel  set  Apparatus 

The  development  of  validated  derailment  criteria  has  been  inhibited 
previously  by  practical  limitations  in  obtaining  a  comprehensive  set  of 
experimental  measurements  of  critical  variables  under  a  range  of  conditions. 
The  apparatus  used  in  this  research  is  designed  to  provide  detailed  measure¬ 
ments  of  wheel  set  motions  and  loading  during  the  evolution  of  derailment 
events.  It  is  also  designed  to  provide  sufficient  control  of  test  conditions 
so  that  the  reproducibility  of  results  may  be  established  and  the  complete 
spectrum  of  phenomena  necessary  for  derailment  prediction  be  explored. 

The  apparatus  used  to  study  wheelclimb  derailment  is  a  one-fifth  scale 
model  wheelset  described  in  detail  in  {29].  The  wheels  and  rails  are 
machined  to  unworn  profiles  from  a  polycarbonate  resin  material  which 
assures  proper  scaling  of  the  wheel/rail  contact  (creepage)  forces  in 
relation  to  the  applied  axle  loading.  The  wheelset  has  lateral,  vertical, 
yaw,  and  roll  degrees  of  freedom,- and  is  attached  to  a  force  measuring 
system  with  vertical  and  lateral  freedoms,  the  latter  simulating  a  generalized 
truck  mass.  The  track  structure  is  tangent  track,  rigid  in  bending  and 
torsion  to  eliminate  interactions  between  wheelset  and  track  deflections, 
and  maintained  to  Class  6+  geometry  [31], 

The  apparatus  consists  of  the  following  elements,  shown  in  Figure  B.l: 
a)  Wheelset  or  truck  model,  dynamically  scaled  from  the  full-scale 
prototype.  Applied  static  and  dynamic  forces  are  scaled  so  that  dynamic 
behavior  of  the  prototype  is  reproduced. 
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Figure  B.l 


b)  Linkage  and  gimbal  systems  providing  the  model  with  vertical  and 
lateral  displacement,  yaw,  roll,  and  pitch  angular  degrees-of-freedom.  The 
sequence  of  freedoms  provided  is  such  that  all  possible  body  orientations 
and  trajectories  are  obtainable,  assuring  normal  wheel  contact  at  all  times. 

c)  Idler  carriage  traveling  along  the  track  with  an  independent 
suspension  system  to  provide  an  attachment  point  for  the  linkage  system. 
Running  on  well-aligned  smooth  rails,  this  rubber-tired  carriage  provides 
a  reference  for  all  force  and  displacement  measurements.  It  also  provides 
a  platform  for  application  of  suspension  forces,  and  is  a  mechanical  link 
to  external  propulsion,  power,  data  reduction,  and  data  recording  equipment. 

d)  Powered  general  purpose  test  carriage,  providing  velocity-controlled 
propulsion,  instrumentation  power,  signal  conditioning,  and  data  recording. 
The  powered  carriage  runs  on  an  independent  rail  system. 

In  the  following  paragraphs  the  wheel  set  and  support  system  are 
described  in  detail. 

Wheelset  Model  and  Gimbal  Syet&n 

The  wheelset  consists  of  two  polycarbonate  resin  wheels  machined 
by  numerically-controlled  lathe  to  the  specified  profile  (Figure  B.2), 
mounted  on  an  instrumented  axle.  The  gimbal  system,  shown  in  Figure  B.3 
and  B.4,  provides  yaw,  roll,  and  pitch  angular  degrees-of-freedom  to  the 
wheelset  or  truck  body.  Each  rotation  is  supported  by  precision  bearings 
mounted  in  close  tolerance  housings.  The  gimbal  unit  is  designed  to 
minimize  deflections  otheT  than  the  desired  rotations,  yet  be  lightweight 
to  reduce  influence  of  gimbal  ineTtia  on  wheelset  dynamics. 


Figure  B.2  Scaled  Wheels  Machined  to  Profile  of  New  Full-Scale  Wheels. 
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ion  View  of  Wheelsct  Mounted  in  Gimbal  System. 


PARALLELOGRAM 

LINKAGE 


Figure  B.4  Side  View  of  Assembled  Wheelset,  Gimbal,  and  Linkage  System 
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The  base  end  of  a  six-component  strain  gauge  balance  is  fixed  to  the 
linkage/balance  adapter,  constrained  to  move  parallel  to  the  track  without 
roll,  yaw,  or  pitch.  By  so  doing,  any  forces  or  moments  transmitted  through 
the  gauge  to  elements  clamped  to  the  metric  end  of  the  gauge  are  sensed  in 
rail  coordinates.  The  yaw/balance  adapter  is  rigidly  clamped  around  the 
gauge, with  its  other  cylindrical  axis  oriented  vertically. 

The  yaw/r oil  body  adds  the  first  rotational  freedom.  Its  tee-section 
encloses  that  of  the  yaw/balance  adapter;  the  yaw  axis  's  supported  by 
bearings,  while  the  cylindrical  section  along  the  gauge  allows  clearance  for 
♦3  degrees  yaw  rotation.  The  lower  cylinder  of  the  yaw/roll  body  provides 
the  axis  of  rotation  for  roll.  On  the  top  of  the  yaw/roll  body  is  a  yaw 
plate,  parallel  to  a  similar  plate  on  the  yaw/balance  adapter.  These 
plates  are  used  for  establishing  the  nature  of  the  yaw  freedom;  at  one 
end  of  the  plates  are  locations  for  yaw  springs,  while  at  the  other  end 
are  locking  clamps  to  prevent  yaw  displacement  as  dictated  by  the  experiment . 

The  roll/axle  body  allows  the  second  rotational  freedom.  The  lower 
bearing  supports  surround  the  cylindrical  housing  of  the  yaw/roll  body, 
with  bearings  used  for  support.  The  yaw  bearings  pass  through  the  center 
of  the  roll/axle  body  with  clearance.  The  upper  bearing  supports  are  for 
the  axle  rotation  freedom.  In  truck  experiments,  the  axle  is  replaced 
by  the  truck  center  bolster. 

The  assembled  gimbal  system  is  shown  in  Figures  B.S  and  B.6.  The 
strain  gauge  axis,  coincident  with  the  roll  axis,  is  parallel  to  the  track 
between  the  wheel/rail  contact  points  (when  wheelset  is  centered).  The 
yaw  rotation  axis  passes  through  the  center  of  the  wheelset. 


Figure  B.5  Rear  View  of  Wheelset  and  Gimbal  Assembly. 


Figure  B.6  Wheelset  and  Ginbal,  Showing  Yaw  Protractor  and  Locking  Device. 
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Lateral /Vertical  Linkage  System 

The  lateral /vertical  linkage  system,  shown  in  Figures  B.4  and  B.7 
connects  the  strain  gauge  balance  to  the  idler  carriage  bulkhead  in  a  manner 
that  sustains  its  orientation  parallel  to  the  rails.  The  bulkhead  is 
rigidly  fixed  to  the  idler  carriage,  so  it  is  a  fixed  inertial  plane  normal 
to  the  rail  x  axis.  Lateral  and  vertical  freedoms  with  no  yaw  or  pitch 
rotation  are  provided  by  the  parallelogram  linkages.  The  linkages  by 
themselves  do  not  provide  adequate  roll  stiffness  to  the  gauge  balance  base; 
this  stiffness  is  provided  by  a  torque  tube  connecting  two  automotive 
universal  joints. 

Idler  and  Powered  Test  Carriages 

The  track-mounted  idler  carriage  with  installed  linkage/gimbal /wheel set 
assembly  is  shown  in  Figure  B. 8.  The  carriage  serves  three  functions: 

a)  Provides  a  fixed  reference  for  all  force  and  displacement  measure¬ 
ments. 

b)  Provides  a  platform  for  application  of  static  and  dynamic  suspension 
vertical  and  lateral  forces  and  moments  to  wheelsets  and  trucks. 

c)  Provides  a  link  to  main  overhead  carriage  for  propulsion,  on-board 
power,  data  reduction,  and  data  recording. 

The  carriage  frame  is  fabricated  of  two-inch  welded  aluminum  tubing 
for  high  strength  to  weight  ratio.  At  each  end  of  the  carriage  are  three 
sets  of  pre-loaded  opposing  wheels  that  provide  the  carriage  with  its 
own  independent  suspension.  Each  wheel  is  aluminum  with  specially  molded 
rubber  tires.  The  rubber  on  the  top  and  side  wheels  is  hard  to  provide  a 
stiff  suspension  as  they  roll  on  the  smooth,  aligned  surfaces  of  the  rail 
top  and  rail  mounting  bracket  side  surfaces.  The  rubber  on  the  wheels  riding 
on  the  rougher  underside  of  the  track  support  structure  is  softer.  The 
pre-load  is  adjusted  through  cams  built  into  offset  wheel  axles. 
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Figure  B.7  Parallelogram  Linkage  and  Torque  Tube  Preserving  Rail 
Coordinate  System  Orientation  of  Gauge  Balance  Through 
Connections  to  Bulkhead. 
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Figure  B.8  Wheel  set  Apparatus  Mounted  in  Idler  Carriage. 
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The  test  carriage  is  powered  by  a  variable  displacement  hydraulic  pump 
with  velocity  feedback  to  provide  precise  velocity  control  (speed  constant 
to  1%).  The  carriage  maximum  velocity  is  12.2  m/sec.  The  complement  of 
signal  conditioning  and  recording  equipment  on-board  is  described  in  the 
following  instrumentation  section. 

8*2  Instrumentation 

The  transducers  employed  in  the  wheel  set  experiments  are  listed  in 
Table  B.l.  The  principal  force  measurement  device  is  a  six-component  internal 
strain  gauge  balance,  shown  in  Figure  B.9.  This  device  measures  all  force 
and  moment  components  transmitted  from  the  base  end  to  the  gimbal -mounted 
model.  Since  the  balance  is  maintained  parallel  to  the  track  the  forces 
and  moments  are  measured  in  rail  coordinates.  When  the  wheelset  is  centered, 
the  axis  of  the  gauge  is  colinear  with  the  contact  points,  as  shown  in 
Figure  B.10. 

Translational  and  angular  displacements  are  measured  by  DCDT-type* 
linear  transducers  and  geared  potentiometers,  positioned  as  shown  in 
Figure  B.10.  The  linear  transducers  measure  track  gauge  and  the  lateral 
position  of  the  balance,  while  the  pots  measure  yaw  and  roll  angles. 
Accelerometers  measure  wheelset  and  track  lateral  accelerations. 

Axle  speed  is  not  constant,  due  to  flanging,  so  that  this  variable 
is  measured  by  DC  tachometer  (Figure  B.10).  This  signal  will  also 
indicate  abnormal  test  conditions  such  as  gross  wheel  slippage. 

*  A  Direct  Current  Displacement  Transducer  operates  on  the  linear  variable 
displacement  transducer  principle.  An  oscillator,  rectifier,  and  low- 
pass  filter  built  into  the  transducer  housing  provide  a  DC  voltage 
output  proportional  to  displacement. 
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TABLE  B.l  INSTRUMENTATION  USED  IN  SINGLE  WHEELSET  EXPERIMENTS 


Transducer 


Variable  Measured  Symbol 


Six-Component  Strain  Gauge  Balance 


DC  Generator  (tachometer) 

Geared  Potentiometers  (2) 

DCDT  Displacement  Transducers  (2) 


Servo  Accelerometers  (2) 


Axial  Force 
Lateral  Force 

Vertical  Force 

Yaw  Moment 

Roll  Moment 

Wheelset  Rolling  Velocity 
Yaw  Angle 
Roll  Angle 

Lateral  displacement  of 
wheelset,  relative  to 
track 

Track  Gauge 

Wheelset  lateral 
acceleration 

Track  lateral  acceleration 
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B.3  Data  Acquisition  and  Reduction 

The  system  used  for  data  acquisition  and  reduction  is  shown  schematically 
in  FigureB.il.  All  data  from  the  wheelset  transducers  is  recorded  on  a 
four  track  frequency  nodulated  (FM)  cassette  tape  recorder.  Eight  channels 
are  nultiplexed  onto  each  of  two  tracks,  recording  lower  bandwidth  signals. 

The  higher  bandwidth  accelerometer  signals  are  recorded  directly  on  the 
remaining  two  tracks.  Lowpass  (25  Hz)  analog  filters  are  used  between  the 
transducers  and  multiplexer  (100  samples/s)  to  minimize  aliasing  effects.  One 
channel  of  the  multiplexer  is  dedicated  to  recording  an  identification  signal 
used  to  mark  the  locations  of  calibration  data  and  running  time  records  on  the 
cassette  tape. 

After  a  series  of  test  runs  have  been  completed,  the  cassettes  are 
played  back  for  digitization  of  the  signals.  The  multiplexed  tracks  are 
demultiplexed,  and  all  signals  are  again  lowpass  filtered  pTior  to  digitization 
to  minimizing  aliasing.  The  sampling  rate  is  200  Hz.  Data  is  stored  on 
magnetic  tape  and  reduced  on  an  HP-1000  minicomputer  system. 


Fi pure  R.  11  Wheelset  Experiment  fata  Acquisition  System 
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B.4  Track  and  Support  Structure 
Track  Structure 

To  perform  repeatable,  systematic  experiments  on  the  dynamics  of  rail 
vehicles  requires  careful  control  and  detailed  knowledge  of  the  track. 

Track  irregularities  provide  the  principal  disturbance  input  to  the  truck, 
and  coupling  between  vehicle  and  track  support  structures  may  be  significant. 
Track  in  service  is  subjected  to  intermittent  loadings,  wear,  extremes  of 
temperature,  and  variations  in  other  ambient  conditions  that  lead  to 
irregularities  in  track  geometry,  railhead  profile,  and  rail  surface 
condition  that  may  vary  with  time  and  track  location. 

The  objective  of  the  track  design  in  this  program  is  to  minimize  the 
influence  of  these  effects  in  order  to  isolate  the  important  wheelset 
dynamic  properties.  It  is  emphasized  that  non-ideal  track  characteristics 
are  extremely  important  to  derailment  in  the  field;  in  future  experiments 
the  importance  of  such  effects  may  be  explored  by  systematically  introducing 
them  into  the  model  experiments. 

Consistent  test  track  conditions  are  maintained  through  the  following: 

a)  The  track  structure  is  designed  to  provide  sufficient  static  and 
dynamic  stiffness  so  that  changes  in  track  geometry  under  load 
are  negligible.  Such  high  stiffness  isolates  the  truck  dynamics; 
in  subsequent  programs  track  flexibility  may  be  incorporated  in 
the  experiments  to  quantify  the  importance  of  this  effect . 
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b)  Track  geometry  is  specified  using  statistical  models,  calibrated 
with  survey  data  taken  at  periodic  time  intervals.  Adjustments 
are  provided  for  vertical  and  lateral  alignment,  gauge,  and  cross 
elevation. 

c)  The  railhead  profile  is  machined  to  be  invariant  over  the  track 
length.  Profiles  from  different  track  locations  may  be  sampled 
over  time  to  check  wear  through  use  of  optical  comparators  and 
contact  impression  tape. 

d)  Uniform  surface  condition  of  the  rails  is  maintained  by  its 
indoor  location  and  periodic  cleaning  with  methyl  alcohol  as 
standard  test  procedure. 

e)  Irregularities  induced  by  temperature  variations  (range  of  27°C 
yearly)  are  minimized  by  design;  measurable  effects  are  eliminated 
through  realignment. 

These  design  goals  are  quantified  as  specifications  in  Table  B.2. 

The  remainder  of  this  section  describes  the  design  and  measured  characteristics 
of  the  track  and  support  structure. 

The  support  structure  shown  in  Figures  B.12  end  B.13  is  an  elevated 
dual  channel  beam  configuration  with  cross  webs  for  increased  torsional 
rigidity.  The  structure  is  supported  on  1/2  in.  (1.27  cm)  threaded  studs, 
permitting  vertical  alignment  and  cross  elevation  adjustment.  The  support 
studs  are  welded  to  steel  cross  ties,  which  are  secured  to  the  concrete  floor  with 
thunderstud  bolts.  Rail  lateral  alignment  and  gauge  adjustment  are  achieved 
through  movement  of  rail  mounting  brackets  along  a  steel  plate  bolted  to 
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TABLE  B.  2  DESIGN  GOALS  FOR  SCALE  MODEL  TRACK 


Parameter 


Goal 


Track  Stiffness: 

-  Vertical 

-  Lateral 


0.12S  m  deflection  foT  44S  N  vertical 
load  on  both  rails 


0.125  wm  deflection 
load  on  one  rail 


for  111  N  lateral 


Alignment : 

-  Vertical  6  Lateral  Spectral  density,  in  scale,  as  good 

as  or  better  than  Class  6  track. 

-  Gauge  28?  *_  0.25  mm 

♦  0.25  mm  (0.08  percent  slope) 


-  Cross-elevation 


Cross  Section  of  Adjustable  Track 
Support  Structure  (full  size). 


Figure  B.13  Installed  213  m  Scale  Model  Test  Track. 
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the  top  of  the  channel  beans.  These  adjustment  mechanisms  are  required 
due  to  inevitable  concrete  floor  settlement,  temperature  gradients,  and 
relaxation  of  internal  stresses  in  large  section  members  that  will  tend 
to  distort  and  misalign  the  track.  Measured  stiffness  properties  of  the 
track  structure  are  given  in  [29]. 

Measurements  of  Model  Track  Geometry 

Track  geometry  was  determined  using  statistical  models  in  [29]. 

The  resulting  spectral  densities  for  vertical  and  lateral  alignment,  cross¬ 
elevation,  and  gauge  are  given  in  Figure  B.14  through  B.17  plotted  with 
scaled  track  geometry  data.  [31] 

Railhead  Profile 

The  rail  profile  used  in  the  model  track  is  a  geometrically  scaled 
reproduction  of  a  new  133  lb  CF  and  I  rail.  The  results  presented  in 
this  report  are  specific  to  the  wheel  and  rail  profiles  used,  and  may 
differ  for  other  geometries.  The  1/40  rail  cant  is  built  into  the  model 
rail  profile,  so  that  the  gauge  edge  of  the  railhead  is  vertical. 

The  rail  section  is  shown  in  Figure  B.18.  The  rails  are  machined 
from  annealed  polycarbonate  resin  bar  stock  using  special  cutting  tools 
with  the  specified  profile.  Tool  wear  was  checked  to  monitor  profile 
uniformity  over  the  rail  sections.  Sections  were  numbered  in  sequence 
of  machining;  they  were  installed  in  the  same  sequence  to  ensure  gradual 
profile  changes  (if  any)  over  the  track  length.  Two  contact  profiles, 
identified  as  A  and  B,  were  machined  so  that  a  second  surface  is  available 
at  no  extra  cost  should  the  first  surface  be  changed  by  wear  or  impact 
damage.  Two  V-grooves  were  machined  on  the  sides  as  seats  for  mounting 
spring  brackets;  a  small  extra  groove  is  added  to  identify  the  contact  profile. 
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SPECTRAL  DENSITY ,  IN  /CYCLE/FT 
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Figure  B.14  Spectral  Density  of  Model  and  Scaled  Class  6  Track  Profile 
(Vertical  Alignment)  £313  . 


SPECTRAL  DENSITY ,  IN  /CYCLE /FT 
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SPATIAL  FREQUENCY  ,  CYCLES /FT 


Figure  B.15  Spectral  Density  of  Model  and  Scaled  Class  6  Track 
Alignment  (lateral)  [31], 


L 


J 


SPECTRAL  OENSITY  ,  IN  /CYCLE/FT 
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SPATIAL  FREQUENCY ,  CYCLES/  FT 


Figure  B.16  Spectral  Density  of  Model  and  Scaled  Class  6  Track  Cross 
Level  [31]. 


MODEL 

TRACK 
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The  rail  is  generally  square  in  cross-section,  and  is  considered 
to  be  rigid  as  mounted  on  the  support  structure.  As  described  previously, 
in  the  present  experiments  rail  flexibility  is  not  included,  so  as  to 
isolate  truck  dynamics. 

Preparation  of  Wheel  and  Rail  Contact  Surfacee 

The  establishment  of  experimental  conditions  that  produce  scaled 
wheelset  forces  and  moments  was  found  to  be  critically  sensitive  to  the 
coefficient  of  friction  p  between  the  wheel  and  rail  surfaces.  Considerable 
*  effort  was  required  to  develop  procedures  for  cleaning  and  final  surface 

finishing  to  achieve  uniform,  realistic  adhesion  characteristics. 

Adhesion  limits  p  have  been  measured  to  range  from  0.1  to  0.5  in  field  tests 
using  full-scale  vehicles,  the  limit  being  a  strong  function  of  surface 
condition,  presence  of  contaminants,  and  wetness,  and  a  weak  function  of 
forward  velocity.  It  is  desirable  to  design  experiments  to  operate  at 
the  upper  end  of  the  range,  since  the  effects  of  all  creep  forces  (which 
are  proportional  to  p)  are  more  readily  measured.  The  polycarbonate  resin 
has  a  friction  coefficient  of  0.52  (from  product  literature)  to  0.3  (from 
roller  rig  measurements).  Under  ideal  conditions,  therefore,  it  is 
suitable  for  use  as  a  material  for  wheel/rail  contact  surfaces. 

For  the  initial  wheelset  experiments  the  wheel  profile  had  a  satin 
finish  resulting  from  the  precise  machining  performed  on  a  numerically- 
controlled  lathe.  The  rail  surface,  machined  using  a  specially  contoured 
flycutter,  had  longitudinal  grooves  due  to  imperfections  in  the  cutting 
tool,  and  lateral,  scalloped  markings  due  to  movement  of  th»  rail  past 
the  cutter.  Although  these  machining  marks  were  less  than  25  pm  in  depth. 
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they  significantly  reduced  the  achievable  adhesion  limit.  Even  with 
extensive  chemical  cleaning  to  remove  contaminants,  the  maximum  friction 
coefficient  measured  (in  yaw  moment  tests  in  flange  contact)  was  about 
0.1. 

A  procedure  was  later  developed  to  remove  these  tool  marks  using 
emery  paper  (grit  Nos.  180  and  400)  and  micromesh  diamond  cloth  (no.  1000) 
successively,  each  used  in  machine  sanders  held  in  jigs  to  preserve  the 
original  rail  profile.  The  contact  surface  is  inked  with  felt  markers 
to  highlight  the  machine  marks;  the  surface  is  abraded  until  all  tool 
marks  are  removed  and  a  smooth  finish  is  achieved.  The  wheels  and  track 
are  chemically  cleaned  with  methyl  alcohol  prior  to  each  run,  with  a 
strong  industrial  solvent  used  periodically  to  remove  grease  accumulated 
after  the  track  has  been  unused  for  extended  periods. 

B.5  Description  of  Wheelclimb  Experiments 

Experiments  are  conducted  at  constant  forward  velocity  and  axle  load, 
preset  at  the  beginning  of  each  test.  Lateral  forces  arc  applied  to  the 
model  as  shown  in  Figures.  19  .  Wheel/rail  reaction  forces  are  measured 
by  the  strain  gauge  balance  which  produces  signals  proportional  to  the 
total  wheel/rail  contact  forces  acting  on  the  wheelset  (axle  L/V  ratio). 
Comparisons  with  theory  are  made  using  the  axle  L/V  ratio,  since  separate 
wheel  force  measurements  were  not  made.  Quasisteady  derailment  (Process  A) 
is  measured  by  progressive  lateral  loading  of  the  model  with  a  ball  screw 
drive  motor,  with  vertical  load,  roll  moment  and  forward  velocity  held 
constant.  During  the  dynamic  experiments  (Processes  B  and  C),  a  pneumatic 


B-32 


Figure  B. 19  Application  of  Forces  to  Kheelset  Model. 
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servo  actuated  photoelectrically  at  specific  locations  along  the  200m  test 
track  applies  lateral  forces  to  the  wheelset  model.  The  magnitude  and 
,  duration  of  the  applied  lateral  forces  are  varied  during  the  series  of 

experiments,  as  well  as  vertical  axle  load  and  forward  velocity.  To 
provide  uniform  initial  conditions  for  the  experiments,  the  wheelset  is 
held  against  stops  near  the  opposing  rail  prior  to  application  of  the 
lateral  force.  Responses  are  recorded  on  an  FM  tape  recorder,  and 
subsequently  digitized  for  detailed  analysis  off-line. 

I  Experiments  on  Process  A  were  conducted  at  yaw  angles  over  the 

range  of  -3°  to  +3°,  for  roll  moments  covering  the  range  possible 
without  wheel lift  at  either  wheel.  To  simulate  Process  B,  a  series 
of  experiments  was  conducted  at  constant  yaw  angles  of  3°  to  -39 
by  locking  this  degree  of  freedom.  Experiments  for  Process  C  were 
conducted  using  a  yaw  spring  suspension  with  low  stiffness.  Tests 
were  conducted  on  the  wheelset  from  below  critical  speed  to  beyond  the 
onset  of  hunting.  Lateral  force  pulse  inputs  were  superimposed  on  the 
track  inputs  to  simulate  dynamic  curving. 

Experiments  were  performed  to  determine  the  reproducibility  of  the 
derailment  events.  This  information  is  useful  in  determining  the  extent 
to  which  derailment  processes  may  be  described  by  differential  equations 
in  a  deterministic  manner,  without  including  random  variations  in  parameters 
or  forcing  functions.  The  results  of  these  experiments  were  highly 
reproducible  at  specific  track  locations,  but  exhibited  variation  in 
results  when  the  same  forces  were  applied  to  the  wheelset  at  different 
sites  along  the  track  as  shown  in  Figure  B.20.  Since  the  rails  were 
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Figure  B.  20  Effect  of  Track  Location  on  Dynamic  Wheelclimb  Results 
for  Equal  Lateral  Force  Inputs. 
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machined  to  a  profile  that  is  relatively  uniform  along  the  track  and 
the  rails  cleaned  with  methyl  alcohol  prior  to  each  test,  the  occurence 
of  derailment  appears  to  be  sensitive  to  minor  local  variations  in 
friction  coefficient  and  railhead  profile. 
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APPENDIX  C 

DESCRIPTION  OF  PHASE  TWO  EXPERIMENTS 

C . 1  Purpose 

The  experiments  conducted  during  Phase  Two  of  the  research  program  were 
designed  to  help  resolve  the  following  issues: 

1.  The  time  durations  of  large  lateral  forces  during  derailment  events 
that  were  measured  in  the  Phase  One  experiments  were  significantly 
larger  than  those  predicted  from  simulation.  The  bandwidth  of  the 
instrumentation  used  in  Phase  One  was  limited  to  25  Hz,  which  would 
filter  lateral  force  impulse-like  signals.  The  bandwidth  of  instrumen¬ 
tation  used  in  Phase  Two  was  extended  to  100  Hz,  which  was  determined 
to  be  sufficient  to  pass  the  predicted  lateral  force  signals. 

2.  In  the  Phase  One  experiments,  the  strain  gauge  balance  was  the  only 
force  measuring  device  employed.  This  created  an  ambiguity  in  re¬ 
solving  dynamic  forces  acting  on  the  individual  wheels,  inertial  forces, 
and  externally  applied  loads.  Under  quasisteady  derailment  conditions, 
with  inertial  forces  insignificant,  the  remaining  loads  were  in  equilibrium, 
so  that  they  could  be  resolved  analytically.  Under  dynamic  loading 
conditions,  it  is  necessary  to  have  independent  measurements  of  the 

load  signals.  In  the  Phase  Two  experiments,  an  instrumented  wheelset 
with  plate  mounted  strain  gauges  was  used  to  yield  lateral  and  vertical 
force  signals  for  each  wheel.  A  load  cell  was  used  to  monitor  the 
external  lateral  force  input,  with  an  accelerometer  used  to  sense 
inertial  forces  in  the  lateral  direction. 

In  the  remainder  of  this  Appendix  the  hardware  and  experimental  procedures 
used  in  Phase  Two  are  briefly  documented. 
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I  C.2  Instrumented  Wheelset  Design 

An  instrumented  wheelset  was  designed  to  provide  measurements  of  lateral 
and  vertical  forces  acting  on  each  wheel,  subject  to  the  following  design 
(  criteria: 

1.  Maximum  sensitivity  of  output  signal  to  input  load. 

2.  Minimum  sensitivity  of  output  signal  to  other  loads  (cross- 

sensitivity)  . 

3.  Minimum  sensitivity  of  output  signal  to  load  point  location. 

4.  Minimum  sensitivity  of  output  signal  to  thermal  and  centrifugal 

stresses. 

5.  Minimum  ripple  in  output  signal  when  subjected  to  continuous  load. 

To  determine  the  best  choice  of  strain  gauge  bridge  geometry  the  strain 

field  of  the  wheel  plate  surface  was  mapped,  for  various  combinations  of  load 
magnitudes  and  application  points.  Using  the  mapped  fields  as  a  data  base, 
a  number  of  bridge  designs  proposed  in  the  literature  were  evaluated  using 
computer  simulations. 

For  the  lateral  bridge,  a  constant  output,  two-sided  configuration  pro¬ 
posed  by  ASEA/SJ  was  selected.  This  configuration,  shown  in  Figure  C.l,  is 
described  in  [32J .  The  principal  advantages  of  this  design  were  moderately  high 
sensitivity  and  no  need  for  processing  of  the  signal  downstream  of  the  strain  gauge 
bridge  amplifier  [32].  The  principal  disadvantages  cited  in  [32],  sensitivity 
■  to  thermal  and  centrifugal  stresses,  proved  to  be  insignificant  problems  for 

the  Lexan  scale  model  wheels  used. 

For  the  vertical  bridge,  an  "A+B"  triangular  output  configuration  proposed  by 
(  ASEA/SJ  was  selected.  This  configuration,  shown  in  Figure  C.2,  is  also  described 

in  [32],  Two  bridges  of  the  geometry  shown  in  the  figure  are  required,  mounted 
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45  deg  out  of  phase  with  each  other.  The  outputs  of  the  two  strain  gauge  bridge 
amplifiers  are  added  together  using  the  following  equation. 

Output  *  Max  (|a| ,  ] B] ,  K(|a|  ♦  |b|))  (C.l) 

where  K  is  a  constant  selected  to  minimize  ripple.  In  the  experiments  reported 
here,  bridge  signals  A  and  B  were  recorded  directed,  with  the  operations  specified 
by  Eq.  C.l  performed  by  digital  computer  during  post -experiment  data  reduction. 

The  final  wheel  designs  required  application' of  twenty-eight  strain  gauges, 
fourteen  to  each  side  of  the  wheel.  Gauges  selected  were  Micro-Measurements  EA-41- 
062DN-350,  selected  for  their  small  size,  convenience  of  mounting,  and  compati¬ 
bility  with  the  Lexan  wheel  material.  The  strain  gauge  signal  bridges  were 
connected  to  their  external  excitation  and  amplifiers  via  slip  rings  mounted  at 
each  axle  end.  The  completed  assembly  is  shown  in  Figure  C.3. 

C.2  Instrumented  Wheelset  Instrumentation  and  Calibration 

The  strain  gauge  bridges  were  excited  and  sensed  using  Analog  Devices  Model 
2B31  strain  gauge  conditioners.  Three-pole,  Bessel  type  active  lowpass  filters 
were  used,  with  a  cutoff  frequency  of  100  Hz. 

The  completed  wheelset  and  signal  conditioning  subsystem  was  calibrated  by 
mounting  the  wheels  in  a  machining  lathe  fitted  with  a  specially-designed  rolling 
wheel  fixture.  The  fixture  was  equipped  with  a  wheel  machined  to  the  same 
cross-section  profile  as  the  rail  used  in  the  track  experiments.  The  fixture 
provided  variable  vertical  and  lateral  loads,  over  a  range  of  load  application 
points.  Using  the  lathe/fixture  apparatus,  the  wheelset  could  be  calibrated  under 
running  conditions  similar  to  those  encountered  during  derailment  tests. 

The  calibration  experiments  consisted  of  applications  of  vertical  loads  in  the 


range  0  to  34  n  (derailment  tests  at  22  n  nominal  load),  lateral  loads  in  the 
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Figure  C-4  -  Instrumentation  System  for  Phase  Two  Experiments 
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range  0  to  18  n  (maximum  lateral  impulse  load  subsequently  measured  was  90  n) , 
for  contact  point  locations  from  the  flange  to  the  outside  edge  of  the  wheel . 

From  these  calibration  tests,  the  sensitivities  of  the  output  signals  to 
applied  loads  were  determined.  The  sensitivities  were  found  to  be  independent 
of  contact  point  location  and  cross-load  to  within  ♦_  10%,  as  desired. 

C.3  Instrumentation  and  Data  Acquisition  Systems 

The  complete  instrumentation  system  for  the  Phase  Two  experiments  is  shown 
in  Figure  C.4.  The  system  is  similar  to  that  used  in  Phase  One,  with  the  exceptions 
of  the  instrumented  wheelsets,  higher  bandwidth  strain  gauge  amplifiers,  and  load 
cell  for  measurement  of  external  applied  lateral  force.  Data  was  recorded  on  a 
fourteen  track  Honeywell  5600  FM  tape  recorder.  The  availability  of  more  tracks 
eliminated  the  need  for  multiplexing  (as  in  Phase  One),  permitting  higher  bandwidth 
for  all  signals.  Subsequent  to  the  track  experiments,  the  data  was  processed  by 
lowpass  anti-aliasing  filters  (100  Hz  bandwidth),  and  digitized  at  400  Hz  on  a 
Preston  Scientific  GMAD-1  analog-to-digital  converter.  Data  storage  and  reduction 
were  performed  on  a  Hewlett-Packard  21MXE  Series  1000  minicomputer  system  located 
at  the  Princeton  University  Gas  Dynamics  Laboratory. 
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LISTING  OF  DYNAMIC  WHEELCLIMB  SIMULATION 
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D-2 


t 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PBOSIkB  TO  S3  BOUTS  TBBBB-BXBSBSXOBkL  BOTXOB 
OP  k  SIKH  BBBBLSET,  FLOS  until 
BXSFLkCBB JUT  OP  Til  THICK  PBkflB.  MBS  IltklU 
COKTkCT  CBOUTII  kBD  KkLKSB  CUIP  POBCBS 


PBOSBkBBBD  BT  L.  I.  SVBBT  BT  JIB  BkILEkT 
TBCBRXCkL  IBSBkICB  XBSTITDTB,  TOf TO,  JkPkR 
VBBSJOR  1  BIT,  1982 


1 

BHtklL 

1 

1 

“  “  “  4 

T 

1 

KklTkB 

l 

4 - 

1 

■"* 

♦— 

» 

» 

i 

TUCK 

1 

T 

i 

HIT 

l 

l 

< 

•BBPEiT* 

> 

» 

1 

4-- 

OOBBK 

1 

♦ 

I 

♦ 

*— - - - ♦ 

- ->|  6  BOB  I 


l 

BBPTZ 

-4 

1 

1 

1 

1 

CBEBP6 

l 

1 

1 

1 

1 

- >| 

BklKEB 

1 

•4 

■asoooio 

■BS00020 
•HS000  30 
BRS00040 

bbsoooso 

•BS00060 
■BS00070 
BBS00060 
BBS00090 
■BS00100 
■BS00110 
■RS00120 
BHS00130 
■BS00140 
IHS00150 
•H SO 01 60 
BRS00170 
BBS00190 
•US00190 
•HS00200 
■HS00210 
■HS00220 
■HS00230 
■BS002U9 
■HS00250 
VHS00260 
■HS00270 
■BS00280 
BHS00290 
SHS00300 
BRS00310 
■BS00320 
■BS00330 
VHS00340 
BBS00350 
BBS00360 
■  BS00370 
■RS00380 
BHS90390 
BN500409 
■HS00U10 
■BS00420 
BBS  004  30 
BRS00440 
■RSOOuSO 
BH500460 
BBSD0470 
BHS00483 
BRS00490 
■RS00S00 
BBS 00% 10 
BBS00520 
■BS00S30 
BRS00540 
BHS00SS0 


f 


it 


] 


L 


j 


D-3 


P IDLER 


- >| 


C««««««<«««<«««««0»»»»»»»»»»»»»»»»»»» 

c««<«««««««««««««»»»»»»»»»»»»»»»»»»»> 
c  - - — - - - 

c  l  iftza  i 

c  - - - - — « 

c 

BEAL  BIX.BZZO 

D1REIS10R  PX  (21.13.13)  .PI  (21,  13.  13). 1(7)  .P3{2)  ,»I  (7)  .ADA  (AO)  . 

1  P1(2),P2(2)  ,XPtl«(3),rS(7) 

BXTERRAL  POSRTT 
COBBOB  /ERIK/  PX.PT 
CORROI  /TPRiD/  TEL .TOR 
COHHOB  /LORD/  PE, REX, PTO, REZO 
COHROR  /PORCRS/  P1.P2.P3 
CORHOa  /ICC EL/  IDDOI.PSDDOT 
COHHOa  /TRUCK/  PS.PSXT 

UTI  BPSI.DEL.T.B  /1.0E-02.0.0.0.0, 1.0R-03/ 

DATA  TBAX.BRAX  /.  20,  $000/ 

OAT A  I  PLOT  /»/ 

BATA  IRIX  /S.OB-OA/ 

DATA  XPLAC  /3*0/ 

CALL  BBRAIL 
CALL  KALTAB 
CALL  TRACK 
CALL  XBXT(T) 

CALL  PBORH(T) 

DO  AS  Xs 1,2 

AS  XP  (P3  (X)  .LB.  0.0)  XPLA6  (X)  »1 

XP  (XPLAC  (1)  .IQ.  1  .OR.  XPLAC  (2)  .10.  1)  80  TO  30 
CALL  POERXT(T.r.OI) 

T-0.0 

CALL  STORE  (T.T) 

BPT«0 

BRXTB  (0.10A) 

10A  PORRAT  (•  BRIER  RX SPLAT  XBTISTAL') 


■RS00560 
BHS00570 
BBS00540 
BBS00590 
RBS00600 
BBS00610 
BBS00620 
BBS00630 
BBS006AO 
■RS006S0 
BRS006S0 
TRS00670 
WBS00680 
RBS00690 
BRS00700 
BRS00710 
■HS00720 
BBS00730 
BBS007«0 
BBS007SO 
THS00760 
THS00770 
BBS00780 
■RS00790 
RBSOOBOO 
BRS00810 
f BS00920 
BHS00830 
BHSOORUO 
RHS0OR50 
THSOOBfO 
TBS00870 
BH50QR80 
■HS00B90 
BRS00903 
RBS00910 
■BS90929 
■HS90930 
BBS709U0 
RH SO 09 50 
IBS  00960 
BHS00970 
BBS009B0 
•HS00990 
BHS01000 
RHS01010 
BRS01020 
•HS01030 
BHSOIOiiO 
BHS01050 
BRSOIOfO 
1HS01070 
•RS010B0 
BRS01090 
AHS01100 


nnno  nonnn 


D-4 


MAS  (5,*)  I*T 

■me 

101  poutt  (•  nnunoi  ikxbsm 

20  CIU  PBOtl(T) 

90  25  !■  1.2 

25  1?  (93  (X)  .It.  0.0)  XPU0(X)-1 

C  IP  (XPU6  (1)  .It.  1  .010.  XPLtB(2)  .It.  1)  00  TO  2* 
C  Cllt  StOtt  (Till.  T) 

C  CO  TO  30 

26  III— 123*5 

■•tixii(a.ntx) 

Ctu  »*1 1  (PDIIXV.T, T,9T,7, tPSI, I, tOX, lit) 

IP  (Xtt  .IQ.  100)  00  TO  20 
TBB0T-9T(2) 

P5P90T*DT  (5) 

IP  (BOD(IPT.IPIOT)  .IQ.  0)  Cttt  STOII(T.T) 

XP  (T  .OT.  TUX)  SO  TO  30 

XP  (IBS  (T  (1) )  .It.  0.0200)  00  TO  35 

XPL1G  (3)  *  1 

CALL  STOtt  (T.T) 

00  TO  30 
35  BPT-BPT»1 

XP  (OPT  .61.  lUX)  SO  TO  30 

IP  (BOD(BFT.IBT)  .ft.  0)  00  TO  20 

■  tITI  (6.100)  fPT.T.I.  T  (1)  «T  (2)  .  T  (!)  ,T  (6)  .T  (7) 

100  PORR1T  (X»,2E9.2,«X,5B10.2) 

60  TO  20 
30  COSTXfOE 

ItXTt  (6,102)  XPLIG 

102  PORBITC  PUSS  -  *,3X10) 

STOP 

EBD 


I  HIT  I 


SBBtOBTIBE  HIT (T) 

tlM)S  It  SISOUTIOI  COSTIOl  PtURBTIRS  PtOR  TttlIRtt 
SETS  PIOCItR  PtttBZTBIS  ft  OB  BtTI  STtTIHEITS 

BEIL  XTI, XXX, XXX .IS, SO, 10. SXX.OZZO.RT. IPIX.KX, RPSX ,HT 
OiniSXOI  I  (?)  ,P3  (2)  ,llPlt  (2)  ,L  (2)  ,IB0  (2)  ,!T  (2) ,  tTR  (2)  ,BTR  (2)  , 

1  CtX  (2). PI  (2)  ,92  (2). PS  (7) 

CORROI  /RPtltR/  If .ITT.XZX.IXX.RT 

CORIOR  /SPitIR/  IT.BT, BZ,BZ, IPSX, tPSZ,tOSB6t,tPBX,BPBX 

CORROI  /TPtlD/  111,101 

COIROB  /RiTEIL/  00, ZB 

CORROI  /lOkB /  PX.IZX.PTO, IXZO 

COIROB  /OBO/  Z,9Z,B0X,PtX,9PII,90PIX,ttPBt,L,tB0,t1. RTB.tTB.CHX 
COHBOB  /POICIS/  PI ,92, P3 

cohroi  /teeny  tidot.psddot 

COIROB  /TBOCR/  PS.PSIT 
IZ10  (1,*)  at,  BT 


VBS01110 
■  ISO  1120 
BBSO 11X0 
BBSOmo 
BBSO 1150 
BBS01160 
BBSO 1170 
BBS01180 
BBSO 1190 
BBS01200 
BBSO 1210 
BBSO 1220 
BBSO 1230 
BBS0 12*0 
BBS01250 
BBS0 1260 
BBSO 1270 
BHS012R0 
BBS01290 
BH501300 
BBSO 1310 
■HS01320 
BBSO 1330 
RHS01340 
BBSO 1350 
BBSO 1 360 
BBSO 1370 
■HSD1380 
BBSO 1 390 
BBSO 1*00 
BBS0 1*10 
BRS01420 
BBS0 1*  30 
BBSO 1**0 
BBSO 1*50 
BBSO 1*60 
BBSC 1*70 
BBS0 1 «  80 
BBS0 1*  90 
BBSO 1500 
BHS01510 
BBSO 1520 
BRS01530 
BBSO 15*0 
BBS0 1550 
BBSO 1 560 
BRS01570 
BH501580 
BBSO 1590 
BBS01600 
BRS01610 
BBSO 1620 
BBS01630 
BBSO 16*3 
BBS01650 


nnnnnnnnnnn  onnno  non  nonnn 


D-5 


BElO  (1,*)  X1I,XXT,XSZ 
IBID  (t,  •)  II,  BT 
HID  (1,*)  II.U 
IE  ID  (1,*)  nsi.DSI 
■  BID  (1,*)  INl.INI 
IBID  (1,*)  B0BB61 
IttD  (1,*)  10,10,1 
■BID  (t,«|  TIL 
■no  (»,•>  pi,  six 
■110  (1,*)  PT0,RM0 

■no  (i,«)  psit 
SI-«.0«B/(3.  0*  <1. 0-B0**2) ) 
cm  cboi(O.O) 

■0-Bf  (2) 

T  (1)  *0.0 
1(2)  -0.0 
TO*-?BL/B0 
r  (3)  *io* 

T  (*)  -PSIT 
T  (5)  *0.0 
T(6)-0.0 
T(7)-0.0 
YDDOT-0. 0 
IS  (2) -0.0 

cm  sbob  (i(i) ) 

■ET0IR 

■ID 


I  TUCK  I 


SOBRODTIBE  TUCK 

■  BIDS  SIB0L1TBD  01  BE1S0IED  TUCK  SBOHETBt  XIP0TS  PBOR  CHS  PX1ES 

BBTOIR 

BID 


I  rOBBXT  | 

♦ - - — 

SDBROOTXIE  P0EIXT  (T.1,01) 

SDBP00TIBB  C1LLBD  «I  D0EBK  TO  BT1L01TB  PXIST  DEBXT1TXTES  OP  til 
STUB  T1IX1BLBS  IT  B1CB  TXBE  STBP.  BQ01TX0BS  1BE  XI  THE  POPH: 
01  (X)  •  (PIB(X)  ♦PBOL(I)  ♦  PS(I)  ♦  rEXT(I))  /  BISS  (I) 

BBBPE 

PBB  *  BBBEL/11XL  XBTBBICTXOB  P OICIS 
PBOL  -  BOLBB  POBCBS 

PS  -  SOSPBISXOI  POBCBS 
PBXT  -  BXTBBBll  POBCBS 

■1S5  >  BBBBB1LXBBD  XBBBTX1  POB  B1CB  STUB  T1B11BLE 


■IS0 1160 
■BS01670 
■BSD  1680 
BPSO 1690 
•BS01700 
■HS01710 
HHS01720 
IB  S0 17  30 
•HS017U3 
■IS0 1750 
•BS01760 
IRSO 1770 
IRSO 1780 
•RS01790 
BHS01603 
BHS01B10 
■BS01R20 
■HS01830 
■RS01840 
•HS01850 
■HS01R60 
■RS01R70 
IRSO 1880 
IRSO 1890 
■HS01900 
■HS01910 
■RS01920 
IBS  01933 
■RS019I0 
IRSO 1950 
■  BSD  I960 
IRSO 1970 
IHSO 1980 
■HS01990 
TRS02000 
IBS02010 
BRS02020 
IHS02030 
IRS020«0 
■R502050 
■HS02060 
IRS02070 
■HS02090 
IBS02090 
•HS0210? 
■RS021 10 
•HS02120 
BHS02130 
■RS02140 
BRS021S0 
BBS02160 
BHS02170 
■HS02180 
■HS02190 
BBS02200 


uuuuu  uuuuu  uuu 


c 


5 


10 

C 


RBU  MD.XTI.XXZ.IXX.BBSS.IS.ST.BXX.IXXO 

BISIISXOI  7(7)  ,BT  (7)  ,P3  (2)  ,11X15  (2>  ,UnS|2|  ,»TO  (2)  ,  MSI* (2)  , 

1  FI  (2)  ,F2(2)  »CBI  (2)  .Ft  8(7)  ,FB0L(7)  ,FS(7)  .FBXT  (7)  , 

2  II.PIB  (2)  «  BIO (2)  ,17(2)  ,  IM  (2)  •  STB (2)  ,1155(7)  «SXCB  (2) 

COBBOB  /8B0/  S»Bt»BBX»FHX,BFIX»BBPIX»ItFIX«I.»lH0»l7,  1TB,  BTH, CHI 

C0BB0B  /BFBBBH/  BB, ITT. XU. XXX, IT 

COBBOB  /CBTBOl/  BBBCB.STBF.B.I.SPS.FX 

COBBOB  /VFBBO/  1U,«t 

COBBOB  /VOID/  FX.BXI.FTO.BZZO 

COBBOB  /FOBCIS/  F1.F2.F3 

COBBOB  /kCCtt/  TBBOT.FSBBOT 

COBBOB  /MOCK/  FS.FSXT 

C1U  8108  (I (1)  ~TBK) 

CKU  HBBTZ  (BBXXS.BBXXS) 

Cltl  CBBCP6  (7.  Bill  5,  B1XX5,  710, 1I0U.  SIOB) 

C1U  KllKIi<lBXIS,BBXXS,?IO,iBM.t,ClI.T.IXeB> 

Clll  FBBBl(FBR) 

CILL  FB0X.BI  (T.FBOL) 

C1U  FSOSF(T) 


BBSS(1)«1.0 
BISS  (2)-Bl 
BBSS  (3)  “XIX 
BBSS  (4) -1.0 
BBSS  (5)  “1X2 
BBSS  (6)  -1.0 
BBSS  (7)“HT 
BO  5  1-1.7 
FBXT  11) -0.0 
FZXT  (7)  “FT0 
FBXT  |5)  -BZX0 
BO  10  X- 1,7 

BT  (I)  -  (FIB  (X)  *FB0l  (I)  »FS  (I)  ♦FBXT  CD  )  /BBSS  (I) 


Bt  (1)  -7  (2) 
B7  (*)  -7  IS) 
BI  (6) -7  (7) 


BBBDCXB-OBBBB  BODZ1 
BT  (3) -0.0 
BT<B) -0.0 
BT  (S)-O.O 

BBT0BB 

IBB 


nmt 


♦ 

i 

♦ 


S0BR00TIIB  SBBF1B 

SBRPLES  TXBX  SXBULB7X0B  IT  BBXF0IB  TXBB  XBTIB7BLS  SIT  BT  XBXT 

BITOBB 

BSD 


•BS02210 

BBS02220 

BBS02230 

IHS02240 

BBS022S0 

BBS02260 

IBS02270 

BHS02287 

BBS02290 

BHS02300 

IBS02310 

BBS92320 

BRS02330 

1HS02340 

BHS02350 

BRS02360 

BHS02379 

BHS02380 

■HS02390 

BHS02400 

VHSO 2410 

BHS02420 

BHS02430 

■RS02440 

■HS02450 

■BS02469 

BBS02U70 

WHS02H80 

BBS92490 

1HS02507 

BHS02S10 

BBS02520 

■HS92530 

1RS92S40 

■BS92SS7 

■RS92S*>0 

BRS92579 

BRS02580 

f RS92599 

BKS02607 

■FS32610 

BBS02623 

BHS02f  30 

BRS026U0 

■HS02650 

BBS02660 

BHS02670 

BH502R90 

■HS02697 

■HS02700 

BHS02710 

BHS02720 

BHS02737 

OH SO 2740 

18502750 


non  nnn  nnnnn  n—  one*  nnnnn 


D-7 


|  STOBB  | 


SBBB0BT1BE  STOBB  (T,T) 

STOIBS  SkBPLBD  SXBOLtTXOB  UU  XB  CBS  FXLBS  FOB  POTOftE  MOCESSIBG 
BUI  10,  L 

DXBBHSXOH  1(7)  ,13(21  ,UHI  (2)  ,1(2)  .BIO  (2)  ,1T  (2)  ,»TB(2|  ,BTB  (2)  , 

1  CBI  (2)  »P1  (2)  ,F2  (2)  .QP  (2) 

COHOB  /fFBHD/  FBI, V0« 

COBBOB  /IfcTEB 1/  10,  EH 

COBBOB  /GEO/  I,DX.DDX,PaX,DFHX.DDPBI,»l.PBk.X.,HH0.1?.»TH,BTH,aiI 

COBBOB  /FOBCBS/  P1,P2,F3 

BBDXSP*T(1)*1000. 

OBBGfc*!  (3)  /BOB 

Ti  •»!(»)  *1000. 

TtEISP-T  (6)  *1000. 

DO  10  1-1.2 
C1-C0S (kLPHk (X)) 

S1*SX1  (kLFBk  (X)) 

(  OP  (I)  —  (-1)**I*(H0*P2(I)*C1*S1)/(BB*F2(I)*S1-C1) 

XF  (F3  (X)  .IB.  0.0)  QF  (X)  *0.0 

•BITE  (2,*)  T,  BHDISP.TBDXSP,  OBEGk,  TIB,  QP,F3 

BBTOBB 

BHD 


4............... 

I  NBkXl  I 

4— ♦ 


SOBBOOTIBB  BBBkXL 
BBkD  STOBBO  HBBkILk  DkTk 

DIBBBSIOB  TX  (HI)  ,T0X  (HI)  ,TOOX  (161)  ,TPIX  (16 1)  .TOPEX  (161)  , 

1  TDDFBX  (161)  .TtLPBk  (161,2)  ,TL  (161,2)  ,TBBO  (161 ,2)  ,TRT  ( 16 1,  2)  , 

2  TfcTB  (161, 2), TBTB  (161,2)  ,TCB1  (161,2) 

COBBOB  /TGBOB/  TX, TDX.TDDX.TPHI, TOPHI, TDDFHX, TALPHfc, TL, TBHO.  THT, 

1  TfcTB, TBTB, TCBX 

COBBOB  /CBTIOL/  BfcHGl,STEP,B,B,EFS,FI 

B>  BE  LIST  /HBBL/  TZ.TDZ  ,TDD2,TFBI,  TOPHI,  TODFHI,  TfcLFBfc  ,TL,TBBO,TBT, 
1  TfcTB,  TBTB,  TCBX 

IBID  (3, Bill) 

DO  10  2*1,161 
TUI  (J)*- TPBI(J) 

TO  FBI  (3)  —  TOPHI  (J) 

TODFHI  (J)  — TDDFBX  (J) 

BBTOBB 

IBD 


4— 

|  KfcLTtB  I 


■HS02760 
BBS02770 
•BS02780 
HHS02790 
BBS02B00 
BHS02810 
VBS02B20 
•BS02R30 
•HS028U0 
•HS02BS0 
BBS02R60 
BHS02870 
BBS02SH0 
BBS02S90 
•BS02900 
BBS 029 10 
•BS02920 
•HS02930 
BBS0299 3 
BH  SO  29  SO 
BBS02960 
■HS02970 
•RS02980 
WHS02990 
HHS03000 
■BS03010 
BHS03020 
■BS03030 
BHS03093 
BHS030SO 
BHSO 3063 
HHSP3070 
BBS033B0 
BHSO  30  90 
BHS03100 
BHSO 31 10 
BBS03120 
BBS03130 
BHSOSm 
•HS03150 
•BS03163 
■HS03170 
•BS031B0 
BHS031 90 
•HS03200 
BBS03210 
BHSO  3220 
•BS03233 
BBS03293 
•BSO3250 
BHSO  3260 
•BS03270 
BBS332B0 
BHSO  3290 
BHS0J300 


D-8 


c  * - - — ♦ 

c 

inioinn  iam> 

bibebsiob  n(2i,u.ixi.n(2itii(ii),mur(ii,M),miif(i).i]| 
cobbob  /uu/  n,rt 
c 

c  aaio  ealkbb  nans  raoa  cas  rxus  (rxta  aos.  si  tbboogr  7ij 

c 

BO  10  X*1,21 

bbb*so«i 

bead  (bbh,200)  ((pxtbbp (k,b)  ,b*i,13),e*i,i3), 

i  t(rxTEHP(x,H),a*i,i3),R-i,i3) 

200  POBBAT  (//////.  13  (/.SX. 13  (P«. 3)  >  ,//////.  13  (/,SX ,13  (P9. 3) ) ) 

DO  20  B*1,13 
BO  20  R-1,13 
pi(i,k,hi*pxtebp(k,s) 

20  ri(X,X,S)*PTTBBP(R,H) 

10  COBTXBOB 
axTota 
IB0 


c 

c  * - - 

C  I  OX08  | 

c  - - - 

c 


sobboutxbe  cxoa(t) 


c 

C  COHPOTKS  COBTACT  0I0BZT1T  PAIABXTEIS  OS  PBBCTXOB  OP  LATEBAL 

c  postxob  op  bbbelset.  nti  is  btobxd  xb  naiua  pobb.  bxtb  positive 

C  DISPL1CBBXBT  OITA  XT  TOE  TOP  OP  TBB  TAB1.I- 


C 

C  I 

C  PBI 

C  DPBZ 

C  BDPBI 

C  * 

C  BX 

C  BOX 

C  MI) 

C  ALPHA  (X) 

C  BP  (I) 

C  BHO(I) 

C  ATB(I) 

C  BTB  (I) 

C  CHI  (I) 

c 

c 

c 


BBXELSET  LATIBAL  POSXTXOB 
BBXELSET  BOIL  A  BOLE 
B(PBX)  /  BT 
02  (PHI)  /  012 

BBXELSET  C.H.  TSXTXCAL  BXSPLACEBBRT 
0  (X)  /  DT 
02  (X)  /  OX 

aOBXXOBTAL  BXSTABCE  PBOB  C.H.  TO  COBTACT  POIBTS 
COBTACT  AB6LE 
B0LLXB6  BA 0X05 

BEOBXTBXC  BEAB  OP  BZBTX  COBTACT  ELLIPSE 
A(TBETA) 

B(TBETA) 

BOBBALIXED  EALIEB  BPXB  PABABETEB 
X*1  POX  BXCBT  BBEEl 
1*2  POB  LEPT  BBEEL 


BEAL  L 

BXBBBSIOB  TX  (161)  ,TBX(161|  ,T00X  (161)  ,TPBX(161)  .TOPHI  (161) , 

1  TBBPBX  (161)  .TA1PBA  (161.2)  ,Tl (161,2)  ,TBBO  (161,2)  ,TBT(  16 1,2). 

2  TATB  (161,2)  ,TBTE (161,2)  ,TCBX  (161,2)  .ALPHA  (2)  , L (2)  ,  BBO (2)  , 

3  BT  (2)  ,  ATB  (2)  ,BTB  (2)  ,CHX  (2)  ,ABX  (100)  .TALPB1  (161)  .TALPB2  (161)  , 
•  TL1  (161)  ,TL2  (161)  ,TBB01  (161)  ,TBB02  (161 ), TIPI  (161) .  TBP2  (161) , 
S  TATH1  (161)  ,TATB2(161),TBTB1  (16 1)  ,TBTB2  (161)  ,TCB1 1(161), 


BHS03310 

BBS03320 

BBS03330 

BB  8033*0 

BBS033S0 

BBS03360 

BBS03370 

BBS033B0 

BHS03390 

BBS03400 

BRS03410 

BHS03420 

BBS03430 

BBS03440 

BBS03450 

BBS03460 

BHS03470 

BBS034RD 

BBS03490 

BJ!  SO  3500 

BBS03510 

BBS03520 

BBS03533 

BBS03540 

BB SO  3550 

BHS03560 

BHS03570 

BBS03580 

BBS03590 

BBS03600 

BBS03610 

BH SO  3620 

BBS03630 

BBS03643 

BHS03650 

BHS03660 

BHS03670 

BHS036P0 

BBS03690 

BHS03700 

BBS03710 

BBS03720 

BHS03730 

BBS037U0 

BBS03750 

BBS03760 

BHS03770 

BBS0378C 

BHS03790 

BBS03B00 

BBS03B10 

BHS0  3820 

BHSO3B30 

BBS03840 

BHS03BS3 
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c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 


6  fCBX2(Hf) 

coaaoi  /inoa/  ii,»tlnD:,tnzl(mzlmpBi,TtiHt,n.nia,tiT, 
1  TkTB.TBTB.TCBI 

coaaoa  /aao/  t.Dx.DDZ.pai.DPii.DDpBi.kLPBi.i.iao.iv.kTB.BTB.cBi 
COBBOB  /CBTBOL/  BkBCB, STXP.B.B.BPS 

BQBXVkLXBCX  (TkLPBfc(l)  .TklPBI  (1)  )  ,  (TklPBk(1t2)  ,T11PB2<1)  )  , 

1  crmi,ni(i)),(n(H2),ni(iji, 

2  (T1BO  (1)  .YBB01  (1) »  »  CTBBO  (1i2)  ,T1B02  (1) )  , 

3  ITBV  (1»  ,*t?1  (1»>  .  (TBV  (K2)  .TBV  2f  t|  > . 

«  (TiTB  ( 1)  .T1TB1  (1)  ).,  (TUB  (M2)  ,T»TB2(  1) )  , 

5  (tBTI  (1)  .TBTB1  (1) )  .  (TiTB  (1B2)  ,  TBTB2(  1) ) . 

6  (TCtld)  .YCBX1  (1)  )  .  fTCBI  (162)  ,TCBI2<  1) ) 

I BY EBPOLkTXOB  IS  PXBPOBBED  XT  SBBXOBTXBE  •XT?"  ft OB  IBB  Sl-BkTH 
LIBBkBI .  VBXCB  XS  BiSBZ)  OB  kXTKXB'S  XTXBkTXVE  BBT  BOD 

TTkBLX-T 

B-BkBGE 

*T— STIP 

XV  (T  .CT.  BkBGB)  TTkBlB»BkBGE 

XV  (T  .IT.  -BkBGE)  YYkBLE*-BkBGS 

cut  XVV  (B.ST.TPHI. V, ITkBIX, B.SPS.PBl.kOI.IES) 

CkU  XVP  (t.ST.TDPHI,  B.YYkBLE, B.EPS.DPBI  .kOX.IEB) 

Cklt  XVP  (B.  ST  .TDDPBX  .  B.  TTkBLB.  B.BPS.DDPBX  .BOX  »XEB) 

CkLt  XVP  (I.ST.TZ.  V.  TTkBLB.  B.  EPS.  S.kOX.XBB) 

Cklt  XVP  (t.ST.TDZ  ,  B,  TTkBLB.  B.  BPS.  DZ.kBX.IXR) 

Cklt  XVP  (X.ST.TDDZ.B.TTkBLk.R.BPS.BDE.kOX.XER) 

Cklt  XVP  (B«  ST.  TklPBI.  R.TTkBLE, B, IPS. kll. kDX.IER) 

Cklt  XVP  (B,ST.TkLPR2.a.ITkBLX.R.BPS.kL2.kaX.XEB) 

UPBk  (1)  “”kl1*PBI 
ktPBk(2)-kt2«PBI 

CkU  XVP  (B.ST.TL1 .B, TTkBLB,  B. BPS, 1(1) .kDX.IER) 

Cklt  XVP  (a,ST.TL2.B,TTkBLE,R,BP5,l(2)  .kDX.IER) 

CkU  XVP  (R.ST.TBROI.B. TTkBLB,  B.XPS.BBO(I) .kDX.IER) 

Cklt  XVP  (B.ST.TRB02. B. TTkBLB,  B.BPS.RBO (2)  .kBX.XBR) 

CkU  XVP  (t, ST. TBV1.B, TTkBLB, B, BPS, RV  (1)  , kOI.IBB) 

Cklt  XVP  (B,ST,TBV2,B, TTkBLX,  B.EPS.BV  (2)  ,  kOI.I  BR) 

Cklt  XVP  <B,3Y,TkTB1.B,TTkBLB.R.BPS,kTR(1).iax.XEB) 

CkU  XVP  (B,ST,TkTB2,B, TTkBLB, B,BFS,kTR(2),kBX.XKR) 

Cklt  XVP  (B.SY.TBTSt.V.TTkBLX.B.IPS.BTB  (1) .kOX.IEt) 

CkU  XVP  {a,ST,TBTB2,B,TTkBlB,B,XPS.BTB(2)  .kBX.XXB) 

CkU  XVP  (B,3T.TCHX1«R,TTkBLX«B,XP5,CBX  (1) «kBX .XEB) 

CkU  XVP  (B, ST, TCBX2.V, TTkBLB, B, IPS, CBX  (2)  .kDX.IER) 

BBTDBB 

EBD 


|  IXBTZ  | 

SDBIOOTXBE  BXBTZ  (kkXXS.BkXXS) 

COBP  DTE  DXBXBSXOBS  OP  BXBTZ  COBTkCT  X1UPSE 


VBS03660 

VBS03870 

BRS03880 

BBS03890 

■RS03900 

BBS03910 

BBS03920 

■BS03930 

BBS03940 

an  so  39  so 

BBS03960 

B9S03970 

BBS03980 

BBS03990 

taso*ooo 

•Hsoaoio 

BBS04020 
VRS040  30 
BBS04040 
•RS04050 
■8S04069 
VBS04070 
BBS04090 
BBS04090 
BBSOUlOO 
BBS04110 
VBS04120 
BBS04130 
•BS04140 
BBS04150 
BBS  04 160 
BBS04170 
BBS04140 
BBS04190 
•HS04200 
BRS04210 
BBS04220 
BHS042  30 
BBS04240 
BHS04290 
BBS04260 
VHS04270 
VHS04280 
BBS04290 
BBS04300 
DHS04310 
V8S04320 
BBS04330 
BBS04340 
BBS043S0 
VBS04360 
BBS04370 
VHS04360 
tltSO 4390 
VBS04400 


ryrtn  nnn  nnn  nnnnn 


nil  1,10 

01 BEBSXOB  P3(2|  .BBXXS  (2)  ,BBXXS(2)  ,100(2)  .STS  (2). 0TB |2)  ,CBI  (2) , 

1  01  MB  (2)  ,1  (2)  »BT  (2)  »P1  (2)  *P2  (2) 

COBBOI  /BBO/  S,DS,DDt,PBX,DPBX,BDPRX, BIPBB, L.BBO.BT,  BTB,BTB,CB1 
COBBOB  /BBTBB1/  IV, BB 

cobbob  /tpbbd/  txl.tob 

COBBOB  /FOBCBS/  P1.F2.r3 
DO  10  X*1,2 

XF  JF3  (1)  .LI.  0.0)  10  TO  10 
FBB>  (F3 (X)  *BB0  (X)/1B)  ••(1./3.) 

1BXXS  (X)  “4TB  (X)  *PBB 
B1IIS  (I)  «BTB  (X)  *PBB 
10  COBTXBOB 
BBTVBB 
BBD 


|  CBBBF6  | 


SOBBOOTXBB  CBEEF6 (X.iiXXS. BBXXS, TBD, BBGIE, SXGB ) 

COBPOTE  BOB-DXBEBSIOB41  CBEBPBGES  XB  MIBB  FOBB 

am  i.bo 

DZ BEBSXOB  FBO  (2)  .BBCLI (2)  , BBXXS |2)  , BBXXS  (2) , T  ( 7)  ,CEX  (2 )  , 

1  BV  (2)  , BLPBB  (2)  ,1  (2)  ,MOUI  ,BTB  (2)  ,BTB  (2)  .SXGB  (2) 

COBBOB  /GEO/  X. OB,  ODE,  FIX,  0FBX,  DOFBX,  B1FBB.1.FBO.BP.  BTB,  BTH,  CU1 

COBBOB  /BBTEBL/  BO, EB 

COBBOB  /TP BBD/  TSL.TOR 

COBBOB  /CBTBOL/  BBBGE.ST.B, B.EPS, PX 

DEFXBE  STBTE  TBBXBB1ES 

TDOT-I  (2) 

OHEGB'T  (3) 

PSI«T{4) 

PSIDOT»I  (5) 

COBPOTE  B0BDXBI8S20BB1  CBEEPBGES 
DO  10  X»1,2 

GBBB  1>  (TE1-C0S  (FSX)  •0BEG4*BT  (X)  ♦  (-1)  ••X«FSXDOT*L  (X) )  /TEL 
IP  (X  .BQ.  1)  XX*2 
IP  (I  ,E0.  2)  II»1 

T2—  (1  (1)  *1  (2)  )  •COS  (BLPBB  (XX)  )  •DMX*TD0T/3IB  (BBS  (BLFBB  (1 ) ) 

1  +BLPBB (2) ) 

GBBB 2“  (“SIB  (FSX)  «T2/TI1)  /COS  (BIPBB  (X)  ) 

BBOBCOBO  (X)/(BO*SflBT  (BBXXS  (X)  *BBXXS  (X))  ) 

XX>BBS  (GBBB1MB0BC) 

SXGB  (I)>1. 

XF  (6BBB1.1T.0.)  SIGB(I)  — 1. 

BTB*GBBB2»BBOBC 
IP  (Z.BC.1)  KTB*-XTB 


BHS044 10 
BBS04420 
BHS044  30 
BBS04440 
•BS044S0 
BBS04460 
BBS 04470 
SBS04480 
BBS04490 
■BS04500 
BBS04510 
■HS04520 
BBS04530 
■HS04540 
BBS045S0 
■BS04560 
BRS04570 
BHS04580 
BBS04590 
BHS04600 
BH504610 
VRS04620 
BR504630 
■RS04640 
BHS046SO 
BBS04683 
•H5 04670 
BBS04640 
WHS0469D 
■HS04700 
BHS04710 
•HS04720 
BBS04730 
BRS0U740 

BHS047SO 
BRS04760 
■HS04770 
BBS  047  BO 
BBS04790 
BHS04BOO 
BHS04B 1 0 
BHS04B29 
BRS04B30 
BBS94R40 
BHS04O90 
VHS04860 
BH504870 
BRS04880 
BRS0489D 
BBS  04900 
BRS04910 
BHS04920 
BHS04930 
BBS 04 940 
•HS049S0 
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I 


( 


f 


► 


C  CONTEST  TO  POX1B  POBB  OSXBG  Ul Ell'S  DBPXBXTIOIS 
C 

TID  (X)  •  (XX**2*1T1**2)  **0.  S 

XI  <<IX  .IS.  0.)  .110.  (IT*  .61.  0.))  116X1  (X|  >90. 
IP  ((XX  .10.  0.)  .1*0.  (IT*  .IT.  0.))  1*6X1  (I)  —90. 
XT  (XX  .IQ.  0.)  60  TO  200 
1*6X1  (X)  -1T1*  (IT1/XX)  *160.  /IX 
200  COITXIOE 
10  COITXIOE 
•■TO I* 

IRD 


■RS04960 
■BS04970 
IBS04980 
**S0«990 
■ISO 5000 
*■505010 
TBS0502C 
•BS05030 
■BS05040 
IBS05050 
*■505060 


C 

C  4 - — — - - ♦ 

C  I  HIKE*  | 

C  ♦ - — - - 

c 


SOBIOPTXIE  K1XIII  (IlXX5.01XXS.Tin.ll6XB.CBX.  T.  SXCI) 

c 

C  SEXECTS  STOIIO  XIXKII  T11XX  POE  6XTII  I/O  110  SPX*  CBX 
C  PEBPOIBS  2-DIREISI OI1L  T1BXE  XBTIIPOXlTXO* 

C 

DXBEBSXO*  PI  (21.  13.13)  ,PI  (21, 13,  13)  ,  10B1  (21)  .10B2(21)  . 

1  CBX  1  (21)  ,CBX2  (21)  ,11X25  (2)  ,01X25  (2)  ,T*D  (2)  ,1*6XE  (2)  , 

2  PI  (2)  ,P2  (2)  ,CBI  (2)  ,T  (13)  ,Y  (7)  ,5X6*  (2) .  P3  (2) 

COBBO*  /II LK/  PX.PT 

COBNO*  /POICXS/  P1.P2.P3 
XOGIC1X  CO*01,CO*02 

D1T1  10B 1/0. 0.1.  025,3*1. 075, 2*1. 125,  2*1. 175,  *.6, 2*9. 5, 2*  11. 75, 

1  3*13., **13.9/, 

2  10*2/1.025.1.075.3*1. 125,2*1. 175.2**. 5,9. 5. 2*1 1.75, 2*1 3., 

3  3*13. 9,6*20./, 

«  CBI1/3*. 0..  1925,. 255, .0.. 2375,. 0..23.2*. 0. .  1B3,. 0, .  1385.  .0, 

5  .0*5.. 05,-0, .0075,. 1125, .1375/. 

6  CBX2/2* 10... 1925,. 255, 10... 2375, 10... 230, 10.. 10... 183, 10., 

7  . 1385. 10. ,.0*5,. 06, 10. ,.0875,. 1 125, . 1375, 10./ 

D1TI  T/0. 0,0. 1,0.2, 0.35, 0.5, 0.65, 0.8, 1.1, 1.*.1. 7, 2. 0.2. 5. 3.0/ 

C 

C  PEBPOBB  K1LKEB  TIB  It  1RIBPOX1TXO*  POB  X1CB  COBT1CT  POUT,  X«1,2 
C 

00  100  X>1 ,2 

IP  (P3  (X)  .GY.  0.0)  60  TO  30 
PI  (I) >0.0 
P2  (I)-O.O 
60  TO  100 
30  comm 

10B>11IX5  (X)  /B1XXS  (X) 

DO  10  J>1, 21 

COBS  1*  (1  OB  .  6T.  10B1  (J)  )  .1*0.  (10B  .XX.  10B2  (J) ) 

C0RD2>  (CBX  (X)  . 6T.  CBX  1  ( J) )  .110.  (CBX  (X)  .XB.  CBX2(J)) 

IP  (COID1  .110.  C0B02)  60  TO  20 
10  COITXIOE 
IPUE-0 
60  TO  999 
20  BPUE-J 
C 


■B505070 
■BSOS080 
•BSOS090 
•BS05100 
■BS05110 
■BS05129 
■BS05130 
■BS051 *0 
■BS05150 
■HS05160 
■RS05170 
•BS05180 
■HS05190 
•8S05200 
•BS0521D 
BB505223 
■BS05230 
BBSD524D 
■R SO  5250 
BBS05263 
■HS05270 
BBS05280 
•BS05290 
■BS05300 
■BS35313 
■RS05323 
•BS05333 
•HS05343 
IBS05353 
■HS05363 
■BS05370 
•BS05380 
■BS05390 
•BS05U00 
■B SO 5* 13 
IBS05U20 
BBS354 30 
0BSO5440 
WBS05453 
•BS05460 
BBS 3 54 7 3 
BRS05483 
•BS05490 
BBS0550P 


( 


L 


J 


nnonn  nnnnn  n  «nnn  non—  n  non 
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c  sbtop  »m  to  innnimoa  booties.  p  xs  ntcnnu  ncuien  z* 
c  van  mu.  xz  xs  ni  coobdzbatb  op  iui  pozbt  pob  xbtbapolatxob 

c 

VB-ABIB1  (TBO  (I), 3.  0) 

XI-12 

DO  25  X 1*1*12 

25  ZP  ( (IB  .OB.  T(I1))  .BBS.  (TB  .IT.  f(X1«1)))  IX-I1 
P-  (?B-f  (XZ) )/(f  (XX*1)-T  (XZ) ) 

CBI1  2~DXflBBSIOBAL  IBTBBPOLATXOB  BOOTHE 

CBlt  XBTP2D  (PI,  J, XZ, IBOtX  (Z)  , P,  1 ,  PP1 ) 

CBU  IBTP2D(PT,J,II.ABCL*(I)  ,P,2.P72) 

P1(X)«>SX6B(X)*PP1 

P2(Z)~PP2 

IP  (I.SQ.1)  P2(Z)*-P2(Z) 

DO  COBTIBDE 

BBD  OP  BALREB  TABLE  I BTEBPOLATZOB  LOOP 
BETUBR 

IBOBT  IBD  -  BO  RALRZB  PILE  POOBB 

19  DISP*1000.  •!  (1) 

BBITE  (2.1000)  OISP.  BTSTEP 
■BITE (5, 1000) SISP, BTSTEP 

1000  POBBBT  (#  BO  DATA  PILE  POBBD—PBOCBAB  TZBBXBATED  AT  *.P7.2, 

V  DISPLACBBBBT  (BB) • ,/,T*2, ' ZB  TEE  '.ZA.’TB  TIBB  STEP*) 

BETOBB 

IBD 


|  PBOBB  | 


SOBBOOTIBI  PBOBB  (I) 

CALCULATES  BOBBAL  POBCB  P3 (I)  AT  SICB  COBTACT  POZBT  PBOB 
TBO-DZBEBSIOBAL  POBCE  IBD  BOBBBT  BALA BCE  OB  BBEELSBT,  IRCZ.ODXRG 
DT BARIC  TBBBS 

BEAL  L,BB,XTT,XZZ,XZX,R0,BT,112,KT,RPRX,RZ,RPSZ,BXX,RZZ0 
0 ZBBBSZOB  P3  (2)  ,T(7)  ,P 2(2)  ,*A(2)  ,BB(2)  ,  ALPBA  (2)  ,1(2)  ,BT(2)  , 

1  BBO  (2)  »ATB  (2)  ,BTB  (2)  ,CBZ  (2)  ,51  (2)  ,C1  (2)  ,P1  (2)  ,PS  (7) 

COBROB  /CEO/  Z,DZ, DSZ,PBX,PPBZ,DDPBZ, ALPBA, L,BRO,BV, ATS. BTH.CHI 

COHBOB  /SPABAB/  BT.BT, BZ, BZ.BPSZ, BPSI, BOBEBA.RPBI, BPKI 

COBBOB  /BP ABA B/  a«,ZTX,XZZ,ZZI,BT 

COBROB  /LOAD/  PZ,RZX,PT0,BZZ0 

COBROB  /POBCBS/  P1.P2.P3 

COBROB  /BATBBL/  BO, ZB 

COBROB  /ACCEL/  TDDOT.PSDDOT 

COBROB  /TBOCR/  PS.PSZT 


■BS0S510 
IBS05520 
BBS05530 
BBS055U0 
■BS05550 
BBS05560 
■BS05570 
BHS05580 
BRS05590 
BBS05600 
BHS05610 
BRS05620 
BBS05630 
•BS05640 
BRS05650 
■HS05663 
■HS05670 
BBS0S683 
■BS05693 
■BS05700 
BBS05710 
■BS05720 
■RS05730 
BHS05740 
■HS05753 
■H505760 
■HS05770 
■HS0S780 
■HS0579D 
■HS05B33 
■RS0581 0 
■BS05823 
■BS0SB30 
■  H505<»43 
■RS05950 
■BSD5860 
■HS05870 
BBS05883 
BMS9S890 
BBS35903 
BHS05913 
BBS 05920 
BBS0S930 
■BS05940 
■BS05950 
BRS05960 
■RS0S970 
■BS05980 
■BS05993 
BRS06000 
■RS06010 
BRS06020 
BRS06030 
■RS06040 
BHS0A050 
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EBOOT-BZ*PBBOT»BBZ*T  (2)**2 
MBOOT-B2BX*TBBOT«BB2IX*I  (2)  *»2 
l«n*II*IBBOMt*tBOMm 

B-BZX-XIX*PBOBOT-  (XZZ-XTT)  •»  (5)  •*2*P1I-Z?I*I  (3)  *1  (5) 

i  •ini*ni>mi*Bni*i[2) 

BO  10  X-1,2 

B1  (I) -SIB  (ALPBA  (Z)  ) 

10  Cl  (X)— 00S(AL2BA  (Z) ) 

112-1(1)  *1(2) 

B 1-112*  (Cl  (1)-I0*V2(1)*S1  (1)  ) 

B2-112*  (Cl  (2)  -00*22(2)  *S1  (2)  } 

KBtO-0.5*  (If  (1)  *EB  (2) ) 

23  (1)-  (1  (2)  *A*B*BBAB*  (PS  (2)-lB*rDD0T) )  /B1 
23  (2)  -(1(1)  U*B*UIM  (PS  (2)  -■■•TODOT) )  /B2 
C 

C  CBECK  201  23  (I)  >  0.  PBEPBBTS  BBCATITB  23  (Z)  ABB  GZTZS  BAB  DUG 
C 

BO  25  X-1.2 

25  23  (X)-AIAZ  1(23(1),  0.0) 

C 

BZTOIB 

BID 

C 


c  4---— -------- 4 

C  |  2BBU  | 

C  —————— -4 

c 

SOBBOUTIBB  2BBB1  (2 BP 

C 

BEAL  l.SU 


BIBZBSIOB  23  (2)  ,  Pit  (?)  ,  ALPHA  (2)  ,1  (2)  ,  BB  (2)  ,P1(2),P2(2)  ,BB0(2). 

1  ATB  (2)  ,  BTB  (2)  ,CBI  (2) 

COBBOB  /CEO/  X.DX.DDZ.BBI.DPBI.BDPBI.llPBA.i.BBO.BB. ATB.BTH, CBI 
COBBOB  /BATBBL/  BO. ZB 
COBBOB  /20BCES/  P1.P2.P3 
BO  10  J-1,7 
10  2BB(J)-0.0 

BO  20  Z-1.2 

2BB  (2|  -PBB  (2)  4PJ  (Z)  •  (B0*P2  (I)  *C0S  (A12BA  (Z)  )  «  SZB  (ALPBA(Z)  ) ) 

2BB  (3)  -20*  (3)  -B0*P1  (Z)  *23  (Z)  *1B 1 1) 

20  21*  (5)  -111  (5)  4  (- 1)  **Z*21  (Z)  *23  (Z)  *1  (Z) 


BETOBB 

EBB 

C 

c  *————— —4 

C  |  PZOLEt  | 

C  4— —————4 

C 


BOBBOOTZBE  2BQ1SB (T.2I01) 

C 

C  COBPOTZS  20BCBS  BBS01TZB6  2B0B  BOLBB  TEBBS  ZB  BQOATZOBS  OP  OOTZOB 
C  ZBV01VZBG  CBOSS  2B0B0CTS  02  AB601AB  BZS2LACEBBBTS  ABB  BE 10 CITIES 
C 

BEAi  i. izz.zrr.zzz .  to,  it,  ft,  ir.ini,  ft , eesj 

BZBEISZOB  1(7)  ,2BB1(7)  .  1(2)  .  *121 A  (2)  ,110  (2)  , B?  (2)  .ATI  (2)  ,BTB  (2), 


BBS06060 

BBS06070 

BBSOfcOBO 

BBS06090 

BBS06100 

BHS06110 

BBS06120 

BBS06130 
BBS061UO 
BBS06150 
BBS06160 
BBS06170 
BBS06180 
BHS06190 
BBS06200 
BBS062 10 
BBS06220 
BB506230 
IBS06240 
1HS062SO 
BBS062S0 
BBS06270 
f HS36280 
BHS06291 
BBS06300 
■RS06310 
BHS06320 
BHS06330 
BRS06340 
BHS063S9 
BBS06360 
B8SP6370 
■BS06360 
BBS06393 
BRS064 00 
BBSOfctt 10 
IBS06U20 
BRSOfcU  30 
BBS084U0 
BBS06450 
BHS06460 
BBS0R470 
BHS06480 
BBS06490 
■R506500 
BHSPfSIO 
■HSOtS20 
BBS06S30 
BBS065U0 
BBS06550 
•BS36560 
BBS06570 
BBS06SB0 
BBS06590 
■RS06603 


UUUOO  uuuoo  uuuuu  uuo  uuuu 


i  cai(2) 

COHO*  /BIO/  X,DZ,0DX,PBI,DPBX,BDPB1,*LPB1,L,BB0,BV,ITB.BTH.CBX 
COBBO*  /SPUIB/  •I,KT,BE,CZ,KPSI,BPSI,B0B£6I,KPBZ,BPflZ 
COHO*  /BPtllH/  BI.ITT.XXX.XXX.BT 

combo*  /iitnv  ao,a* 

COBBO*  /ICCIL/  TDDOT.PSDDOT 

p*dot-dpbx*t(2) 

BO  10  0-1,7 
10  PSOL  (J)*0.  0 

PEOL  (3)-  (7  (5)  *PBDOT*PBI*PSDDOT)  •XTT-T  (5)  *PBDOT*  (XXZ-XXX) 

PI01  (S)-T(S)  •PBI*PSDOT*IXI“ITT*T  (3)  •  PH  DOT 
1  - (XXX-ITT) •PBI*PBD0T*I (5) 

BZTDBR 

HD 

1  PSOSP  I 


SOBBOOTXaa  PSOSP (T) 

COMPOTES  POBCBS  BEStlTXBG  PB08  DXSPLICEBMT  OP  PBXBMT  SOSPEBSIOIS 
BETIEE*  BBBELSBT  MO  BOGIE  PB1BE,  XICLODX *6  UTtlll,  1EBTXCIL,  BOLL, 
**D  TB*  STirriBSS  MO  DIBPXBG,  PLDS  DUG  OB  BBMXBG  POBCE  O*  BILE 

BEIL  L,KX,KPBX,KZ, BSPX 
DIBEISXO*  1(7), PS  (7) 

COBBO*  /SPMBB/  BT,IT,B2,KZ,KPSX,BPSl,BOBECi,KPBX,BPBI 

COBBO*  /TBOCK/  PS,PS1T 

PS(1)-0.0 

PS  (2)  -H* (I  (6)  -1  (1) )  ♦BT»  (T  (7)  -1  (2) ) 

PS(3)— B0BZC1*T(3) 

PS  (*)-0.0 

PS  (5)  — BPSI-T  (5)  ♦IPS!*  (PSXT-T  (•) ) 

PS  (6) -0.0 

PS  (7)  -II*  (T  (1)  -I  (6)  )  TII*  (T  (2)  -1  (7)  ) 

BETOBR 

EWD 


|  BLOC*  BIT!  | 


BLOCK  D1T1 

OSED  TO  XIXTItLXZE  PMIBETBBS  X*  BLOCK  COBBO*  /CBTBOL/ 

COBBO*  /CBTBOL/  BIBGE,STEP,B,B,BPS,PX 

BEIL  BMGE/*0.6*B-03/,  STEP/5. 0BB-0B/,BPS/-0.0S/»  PI/3. 1*159/ 
XBTXGBB  B/161/, B/3/ 

EBB 


SHTICE  BOOTIBBS 


BUS 066 10 
BBS06620 
BBS06630 
BHS06640 
BBS06650 
BBS06660 
■HS06670 
BBS06680 
BBS06690 
BBS06700 
BBS06710 
BBS06720 
BBS06730 
«BS067«0 
BBS067S0 
BBS06760 
BHS06770 
•BS06780 
BBS06790 
BBS06600 
BBS 068 10 
BBS06820 
BHS06830 
BHS06940 
BHS06850 
BBS06860 
BRS06870 
BBS06880 
BBS06890 
WBS06900 
BBS06910 
*HS0f 920 
BH506930 
■HS06940 
VH5069S0 
*HS96°60 
BHS06970 
■BS06980 
BBS0699 0 
■BS07000 
BHS07010 
1KS07020 
B8507030 
BHS07040 
BBS07050 
BHS07060 
BBS07070 
■HS07080 
BBS07090 
■BS07100 
BBS07110 
BHS07120 
BBS071 30 
BHS071UO 
BHS07150 
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MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BU»EAu  OF  STANOAADS  -  «»6S  "  A 


uuouuuo 


i  mm  i 


sbbvxcs  loniu  to  runii  Tao-oxaBisioau  xbtkbpoutxob 

sobioot: it  mm (P.  j, xx, iku. p, xnu, mot) 
BXBBBSIOB  P(21. 13,  13)  ,1011  (IS)  ,1012  (IS)  ,101  (30) 
un  IPS/- 1.  IB-02/ 

®0  10  1-1,13 
toil  (I«1)-P(J,XX,X) 
io  io«2(x*i)-r(j,  (xx«i),x) 

ion  (i)-ro,xx,2)*(~i)**xrue 

1012  (1)-P(J,  (IX *1)  ,2)*(-1)**XPLIG 

1011  (15)-T  (3,XX,12)*(-1)*«XPtlS 

1012  (IS)  -P  (J,  (1X»1)  ,  12)  •  (- 1)  ••IPLIC 

cm  ZIP  (“105.  ,  IS,  ,1011,  15, 11011,3, BPS, PI #101, III) 
cm  XTP  (-105.  ,15.  ,1012,  15,IB6LB,3,EPS.?B. 101,  XBB) 
PliOt- (1.0-P)»Pl  ♦  P*PB 

izroxi 

110 


IRS07160 
■BS07170 
■BS07180 
IRS07190 
•BS07200 
■IS07210 
■BS07220 
■BS07230 
BBS072B0 
BBS07250 
BR  SO  7  260 
BBS07270 
BBS072B0 
BBS07290 
BBS07300 
BBS07310 
BBS07320 
■BS07330 
BBS073B0 
BBS07350 
IBS07360 
BHS07370 
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APPENDIX  E 

PATENT  DISCLOSURE  STATEMENT 

In  accordance  with  the  patents  rights  course  of  the  terns  and 
conditions  of  this  contract,  and  after  comprehensive  review  of  the 
work  performed,  it  was  found  that  no  new  patentable  items  were 
produced  tinder  this  contract.  However,  significant  innovations  and 
advancements  to  the  field  of  research  on  derailment  safety  were  achieved. 

The  following  hardware  and  software  items  were  developed  as  a 
result  of  the  sponsorship  of  research  under  Contract  DOT-TSC-1603: 

1)  Scale  model  wheelset  experimental  apparatus  for  use 
in  studies  of  derailment  processes. 

2}  Computer  programs  for  analysis  of  measured  wheel/ 
rail  forces  for  comparison  with  proposed  derailment 
criteria. 

3]  Measurement  and  data  reduction  methods  for  on-line 
diagnostics  of  impending  derailment  conditions. 

4)  Computer  programs  for  simulation  of  quasisteady  and 
dynamic  wheelclimb  derailment  processes. 

The  above  items  are  potentially  useful  to  both  model  scale  and  full 
scale  experiments  and  analysis  of  rail  vehicle  derailment  safety. 


DATE 
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